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Abstract. 

The LIT method has allowed ab initio calculations of electroweak cross sections 
in light nuclear systems. This review presents a description of the method from both 
a general and a more technical point of view, as well as a summary of the results 
obtained by its application. The remarkable features of the LIT approach, which 
make it particularly efficient in dealing with a general reaction involving continuum 
states, are underlined. Emphasis is given on the results obtained for electroweak cross 
sections of few-nucleon systems. Their implications for the present understanding of 
microscopic nuclear dynamics are discussed. 
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1. Introduction 



i.i. Preliminary remarks 

A very challenging problem in quantum mechanics is the ah initio calculation of the 
cross section for a perturbation induced reaction involving a many-body system. With 
ah initio we intend a calculation that requires a Hamiltonian H and the kinematic 
conditions of the reaction as only inputs, and treats all degrees of freedom of the many- 
body system explicitly (microscopic approach). In general, if the reaction implies a 
state belonging to the continuum spectrum of H, the challenge may become enormous 
since one has to deal with a many-body scattering problem, which may lack a viable 
solution already for a very small number of constituents in the system. One is then 
forced to introduce various approximations that are either based on rather general 
physical considerations or need to be validated by experiment. This situation, however, 
is particularly unsatisfactory in cases where experiments cannot be performed, like 
e.g. for some nuclear reactions of astrophysical relevance, or where the object under 
investigation is just the main input of the calculation, i.e. the potential V in the 
Hamiltonian. This situation is typical for nuclear physics as well as for any 'non- 
fundamental' theory, where one would like to test the reliability of the 'effective' degrees 
of freedom in the Hamiltonian and of the interaction V . In this case the comparison 
between theoretical results and experimental data is expected to give that information, 
but such a comparison risks to be inconclusive if the quality of the approximation is not 
under control. In both cases described above an accurate ah initio calculation may be 
demanded. 

The difficulty in calculating a many-body cross section involving continuum states 
can be understood if one considers that at a given energy the wave function of 
the system may have many different components (channels) corresponding to all its 
partitions into fragments of various sizes. Already in a rather small system of four 
constituents the two-, the three- and the four-body break-up channels contribute at 
energies beyond the so-called four-body break-up threshold. In configuration space 
the task consists in finding the solution of the four-body Schrodinger equation with the 
proper boundary conditions. It is just the implementation of the boundary conditions 
for a continuum wave function which constitutes the main obstacle to the practical 
solution of the problem. In fact, the necessary matching of the wave function to the 
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oscillating asymptotic behaviour (sometimes even difficult to be defined unambiguously) 
is not feasible in practice. In momentum space the situation is as complicated. The 
proper extension of the Lippmann-Schwinger equation to a many-body system has been 
formulated long ago with the Faddeev-Yakubowski equations [HIS]- However, because 
of the involved analytical structure of their kernels and the number of equations itself, 
it would be very hard to solve the problem directly with their help, at energies above 
the four-fragment break-up threshold, even for a number of constituents as small as 
four. The great merit of the LIT method is that it allows to avoid all the complications 
of a continuum calculation, reducing the difficulties to those encountered in a typical 
bound-state problem, where the boundary conditions are well defined and much easier 
to implement. 

1.2. Some historical notes 

The LIT method [3] is the natural extension of an original idea [1] to calculate reaction 
cross sections with the help of integral transforms. This kind of approach is rather 
'unconventional'. It starts from the consideration that the amount of information 
contained in the wave function is redundant with respect to the transition matrix 
elements needed in the cross section. Therefore, one can avoid the difficult task of 
solving the Schrodinger equation. Instead one can concentrate directly on the matrix 
elements. With the help of theorems based on the closure property of the Hamiltonian 
eigenstates, it is proved that these matrix elements (or some combinations of them) can 
be obtained by a calculation of an integral transform with a suitable kernel, and its 
subsequent inversion. The main point is that for some kernels the calculation of the 
transform requires the solution of a Schrodinger-like equation with a source, and that 
its solutions have asymptotic conditions similar to a bound state. In this sense one can 
say that the integral transform method reduces the continuum problem to a much less 
problematic hound- state-like problem. 

The form of the kernel in the integral transform is crucial. The reason is that 
in order to get the quantities of interest the transform needs to be inverted. Since 
it is normally calculated numerically it is affected by inaccuracies, and inverting an 
inaccurate transform is somewhat problematic. Actually, when the inaccuracies in the 
input transform tend to decrease, and a proper regularization is used in the course of 
inversion, the final result approaches the true one for various kernels [5]. However, the 
quality of the result of the inversion may vary substantially according to the form of the 
kernels, even for inaccuracies of similar size in the transforms. In particular, when, for a 
specific kernel, the accuracy of the transform is insufficient, the result may be corrupted 
with oscillations superimposing the true solution. 

In |1] the Stieltjes kernel was proposed and its reliability was tested and discussed in 
simple model studies. Later, in a test of the method on a realistic electromagnetic cross 
section, calculated also in the conventional way for the deuteron [6], it has been found 
that the use of the Stieltjes kernel is not satisfactory, since it leads to quite inaccurate 
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results. The problem with the Stieltjes kernel can be understood if one notices that 
its form is not qualitatively different from that of the Laplace kernel. In fact it is 
well known that the problem of the inversion of a Laplace transform is extremely ill 
posed when the input is numerically noisy and incomplete [7]. Nevertheless the use 
of Laplace transforms is common in various fields of physics, from condensed matter 
to lattice QCD, and elaborated algorithms (e.g. the maximum entropy method [8]) are 
sometimes employed for its inversion. 

The problems encountered in inverting the Stieltjes as well as the Laplace transform 
has led to the conclusion that, differently from those two cases, the 'best' kernel should 
be of a finite range. Its extension should be, roughly speaking, about the range of 
the quantity to be obtained as the result of the inversion. This allows to eliminate 
extraneous low-frequency oscillations in the inversion results, even at a moderate 
accuracy of the input integral transform, while high-frequency oscillations are excluded 
by a regularization (see section [315]) . At the same time, of course, the transform has to 
be calculable in practice. In [3J it has been found that the Lorentzian kernel satisfies 
both requisites and the analogous test, as had been performed in [6] for the Stieltjes 
kernel, has led to very accurate results. 

The present work constitutes the first comprehensive review article of the LIT 
method. Only short summaries of this approach have been published in [9l [TO] . 
Overviews of the LIT method and of its applications have been presented at various 
conferences [HI [12] 

1 . 3. Descriptive plan of the review 

In section [2] it is shown that the integral transform method in general, and the LIT 
in particular, can be formulated for any perturbation induced reaction, of inclusive 
(section 12. ip as well exclusive (section 12.21) character, involving a general system of N 
interacting particles. There it becomes evident that another remarkable feature of the 
method is that the calculation of an inclusive cross section turns out to be simpler than 
the calculation of any single exclusive cross section. Among the exclusive reactions 
we focus in particular on the two-fragment case (section I2.2.ip . since up to now the 
LIT applications to exclusive reactions have been limited to this kind of reactions. The 
application of the LIT method to the case of charge fragments requires some clarification, 
therefore we discuss it in section 12.2.21 The application of the LIT method to other 
reactions like for example those which are not induced by perturbation (strong reactions) 
is discussed in section 12.31 

Once the general idea of the method is outlined, in section [3] we turn to its practical 
implementation. In order to apply the method in an efficient way it is necessary to 
understand better the role of the parameters of the Lorentzian kernels. This is illustrated 
in sections 13.21 and 13.31 after having defined a few key quantities in section 13.11 Some 
particular cases, e.g. energy-dependent transition operators, are discussed in 13.41 The 
inversion of the transform is discussed in section 13. 5[ while in sections 13.61 and 13.71 a 
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more technical description on the practical implementation of the LIT method is given. 

In section H] the LIT approach is compared to other methods that exhibit some 
similarities to the LIT method. 

As already stated above, in all cases the core of the LIT approach resides in a 
non-homogeneous Schrodinger-like equation, whose solution has asymptotic conditions 
similar to a bound state, and therefore can be solved with a bound-state technique. 
In section |5] we describe one of such techniques that has been successfully applied to 
study nuclear, atomic and molecular few-body systems [131 HH [151 [161 [171 HE]. It belongs 
to the class of methods that make use of expansions on complete sets. In this specific 
case the basis set is given by hyperspherical harmonics (HH) functions. The approach 
is called the HH method and comes in two variations, which differ essentially only by 
the way the convergence to the exact solution is 'accelerated'. The first one (known as 
CHH) makes use of a correlation function, which correlates the pure HH basis in order 
to reflect the characteristics of the potential. In the second one (known as EIHH) an 
effective interaction is introduced, which is built from the bare potential by means of a 
similarity transformation. Details on CHH and EIHH are given in Appendices A and 
B. The results presented in section [8] have been obtained using both approaches. 

In section [6] the formalism of the electroweak perturbation induced reactions, 
discussed in section IH are described in greater detail. 

Section [7] contains four numerical tests of the LIT method. In sections 17.11 and 17.21 
two interesting test cases for the reliability of the inversion procedure, as well as of the 
whole LIT approach, are presented. In particular in section [721 the comparison between 
a response obtained via the LIT method and the corresponding LCZR of section [474] is 
discussed. In section [731 a test on the inversion algorithm is performed on ^He and finally, 
in section 17.41 a comparison between the eigenvalue method described in section 13.61 and 
the use of the Lanczos algorithm to calculate the LIT (section 13.71) is presented. 

Section [8] presents a selection of the results of the various LIT applications. Up 
to now the method has been tested and largely applied to reactions induced by a 
perturbation (the electroweak probe) on a strong interacting system (the nucleus). There 
is a reason for this choice, which we explain in the following. 

In section [9] a brief summary of the review is given. 

As already pointed out above, the strength of an ah initio method like the LIT 
resides in allowing to study the effective degrees of freedom in the Hamiltonian and 
the reliability of the potential V . Nowadays this is a much debated topics in nuclear 
physics and the LIT method represents a unique tool to clarify some important issues. 
In particular it can address the question of the nature of the nuclear force. Questions 
like whether the force is exclusively of two-body nature or three- or more body forces 
are required to describe the nuclear phenomenology can be addressed also in reactions 
with more than three nucleons. At the same time one is able to discriminate among 
potentials derived within purely phenomenological or semi-phenomenological (boson 
exchange) or effective field theory approaches. But more than that and differently from 
the purely strong reactions, electroweak reactions with nuclei can provide additional 
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information on the relevant degrees of freedom in nuclear physics. In fact the real or 
virtual photons, as well as neutrinos, probe not only the explicit degrees of freedom in 
the Hamiltonian (protons and neutrons), but also the implicit ones that generate the 
exchange (charge or weak) currents. In the examples reported here particular emphasis 
on these aspects is given. 



2. The LIT method: general theory 

2.1. Inclusive processes 

In inclusive processes the quantities of interest have the following structure 

r{E) = ^ d^{Q\^,){^,\Q')5iE, - E), (2.1) 

where are solutions to the dynamic equation 

{H~E,)\^,) = 0, (2.2) 

and H is the Hamiltonian of the system. The set l"^^) is assumed to be complete and 
orthonormal. 



f/7|^7)(^7l = l- (2-3) 



The integration and summation here and in (12.11) go over all discrete states and 
continuum spectrum states in the set. We suppose that the norms {Q\Q) and {Q'\Q') 
are finite. 

In the case of perturbation-induced reactions one has \Q) = 0|\l/o), \Q') = O'|\l/o), 
where |\l/o) is the initial state in a reaction (generally the ground state of the system 
undergoing the perturbation), and O, O' are transition operators. Then one has 



iE) = ^ rf7(M/o|0^|^-,)(vI/-,|6'|^o)5(i?7 - E). 



(2.4) 



Here is a set of final states. If O = O', the quantity (12.41) may represent a response 
function or, in general, a contribution to the response of the system to a perturbative 
probe transferring energy £^ to it. 

When the energy E and the number of particles in a system increase the direct 
calculation of the quantity r{E) becomes prohibitive (only in case that E becomes 
larger than the interaction energy, perturbation theory can be used). The difficulty is 
related to the fact that in these cases a great number of continuum spectrum states 
contribute to r{E) and the structure of these states is very complicated. 

The approach to overcome this difficulty that is presented here can be considered 
as a generalization of the sum rule approach, since the use of the closure property of the 
Hamiltonian eigenstates plays a fundamental role. Consider for example a simple sum 
rule for the quantity (12. 4p . based on the closure property (12.31) i.e. 



r{E)dE = {Q\Q'). (2.5) 
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The calculation of this quantity is much easier than a direct calculation of r{E) itself, 
since it requires the knowledge of \Q) and \Q') only. This can be obtained with bound- 
state methods since we have supposed that \Q) and \Q') have finite norms. However, 
this sum rule contains only a limited information on r{E). In order to get much more 
information about it we consider instead an integral transform 

$((7) = ^ K{(T, E) r{E) dE (2.6) 

with a smooth kernel K (specified below). This yields 

$(a) = ^ d^{Q\^,)Kia,E,){^,\Q') 

d^{Q\Kia,H)\^,){^,\Q'). (2.7) 



Using the closure property f l2.3p one obtains 

^a) = {Q\k{a,H)\Q'). (2.8) 

Therefore equation (12. 8p may be viewed as a generalized sum rule depending on a 
continuous parameter a. With a proper choice of the kernel K the right-hand side 
of (12.81) may be calculated using bound-state type methods. Once $(0") is available (12.61) 
may be solved to obtain r{E) via an inversion of the transform. 

Our choice of the kernel K{a, E) is such that both the calculation of $(cr) and the 
inversion of (12. 6p are feasible. We choose [3] 

K(a,E) = — 1— -. (2.9) 

^ ' ' {E-a*){E-a) ^ ' 

Notice that the energy parameters a that we consider are complex. For convenience we 
define them as 

a = Eo + aR + iai , (2.10) 

where Eq is the ground-state energy, and o"/ 7^ 0, so that K{a^ E) is actually a Lorentzian 
function centered on Eq + aji, having aj as a halfwidth 

Kian, ar, E) = ^ . (2.11) 

{E- Eq- anY + aj 

Then the integral transform (12.61) becomes 

L(a^, aj) = jfdE /^^^ , . (2.12) 

^ {E - Eo- any + aj 

Here and in the following the integral transform $(o") with a Lorentz kernel is denoted 

by L{aFi, <^i)- Using the definition (12. 4p it is easy to show that the quantity (I2.12p may 

be represented as 

L{aR,ai) = {%^'), (2.13) 
where the 'LIT functions' \E' and \E'' are given by 

1^) = (H-E^-an-iaX' d\^,), (2.14) 



= (H~E,~aR- laA 6'|^o) • (2.15) 
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These functions are solutions to the inhomogeneous equations 

(^H-Eo- an - ^a,) |^) = 0|^o) , (2.16) 

(^H-Eo-ctr- m,) 1^') = O'l^o) • (2.17) 

Let us suppose that 6' = 6 {6 = 6'). In this case L{a) equals to (^|^) ((^'|^')). 
Since for cr/ 7^ the integral in fl2.12p does exist, the norm of |\^) d^E'')) is finite. This 
implies that and are /oca/ued functions. Consequently, (12.161) and f l2.17p can be 
solved with bound-state type methods. Similar to the problem of calculating a bound 
state it is sufficient to impose the only condition that the solutions of f l2.16p and (12.170 
are localized. This means that in contrast to continuum spectrum problems, in order to 
construct a solution, it is not necessary here to reproduce a complicated large distance 
asymptotic behaviour in the coordinate representation or singularity structure in the 
momentum representation. This is a very substantial simplification||| 

Obviously, localized solutions to (12.160 and (I2.17P are unique. Once L{a) is 
calculated r{E) is obtained by inversion of the integral transform with a Lorentzian 
kernel (12.120 ('Lorentz integral transform'). The inversion of the LIT will be discussed 
in section 13. 5[ 

Before ending this section we should mention that if the Hamiltonian is rotationally 
invariant it is useful to expand the states 0\^q) and 0'\^q) over states possessing given 
values J and M of the angular momentum and its projection. Then the whole calculation 
may be done in separate subspaces of states belonging to given J and M. Furthermore, 
the calculations are M independent. An explicit example is discussed in section 16.1.21 



2.2. Exclusive processes 

Here we illustrate the LIT method for perturbation induced exclusive processes. Some 
comments on non-perturbative processes will be listed as well, at the end of this section. 

The relevant matrix element that one has to calculate in the case of an exclusive 
perturbation induced reaction is 

Mf,{E) = {m-^{E)\d\^!,), (2.18) 

where O is a perturbation operator that causes the transition, |\l/o) is again the localized 
initial state, while \^^{E)) represents the continuum final state with energy E. The 
calculation proceeds as follows [ij. We set pS] 

\^-f{E)) = A\ME)) + {E-H- tr])-'\ME)) , (2.19) 

where A\(f)f{E)) is the 'channel' function, A being the antisymmetrizer In general 
this will be the properly antisymmetrized product of fragment bound states times their 

I In the case of Faddeev-Yakubowsky equations one also needs not impose boundary conditions to 
get a solution. However, this is achieved due to an involved structure of these equations that seriously 
complicates calculations already for A = A. In contrast, we use a simpler Schrodinger operator. 
§ Note that here any expression containing x ± 177 is understood to be evaluated at r/ > with the limit 
77 ^ then taken. 
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free relative motion sub-states of given momenta (we temporarily neglect the long-range 
Coulomb interaction between fragments in the final state; the Coulomb case is discussed 
in section [2.2.2p . To deal with states of well defined rotational quantum numbers one can 
write A\(t)f{E)) as properly antisymmetrized product of fragment states and free relative 
motion sub-states of given orbital angular momentum quantum numbers, coupled to 
proper total angular momentum. The states \(f)f{E)) in (12.191) are 

10^) ^ (i7 - E)A\<Pf) ^ AVf\ct>f) , (2.20) 

where Vf is the interaction between particles belonging to different fragments in the 
non-antisymmetrized state 

Using (I2.19p . we rewrite the reaction amplitude (I2.18P as 

Mf,{E) = {A<Pf{E)\d\^>^) + {4>f{E)\{E -H + tv)-'d\<l!o) • (2.21) 

The first term in (12.211) is called the Born term and can be computed directly. The 
second term takes into account the final state interaction and represents the difficult 
part of the problem. However, if we represent it as a sum over the discrete spectrum 
En and integral over the continuous energies E', i.e. in the form 

{ME)\{E -H + tv)-'d\^o) = J2(E - E^)-'Ff,{E, E^) 

n 

POO 

+ / {E- E' + t7])~^Ff,{E,E')dE', (2.22) 

where the form factor Ffi is defined as 

Ff,{E,E') = ^dj{ME)\'^,){^^\d\^o)S{E^ - E') , (2.23) 

we realize that (12.231) has the same formal structure as (12.11) . with E —* E' . Therefore 
the form factor (I2.23P can be calculated at a given E value as a solution to the equation 

$/,(E,a) = ^ K{a,E')Ff,{E,E')dE'. (2.24) 

This equation is similar to (12.61) . Therefore Ffi{E,E') can be obtained calculating first 
^fi{E,a) from (12.81) and then inverting the transform. 

Also in this case we adopt the Lorentz kernel. However, it is useful to express it as 

K(a,E') = — i i ^ , (2.25) 

^ ' ' 2iai \E' -a E' -a* ) ' ^ ^ 



where a = crj? + zcr/ . Then j20| : 

Lf,{E,a) = {2iai)-\^f{E){\i>,) - \i>,)), (2.26) 

where 

\^^) = {H-aR-iai)-^d\^,), (2.27) 

\^>2) = {H -aR + iai)-^d\^Q). (2.28) 
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Therefore and |\^2) are the solutions of the two inhomogeneous equations 

^^-a^-^a,) l^i) = 6|*o), (2.29) 
H~aR + iai'^ 1^2) = 0|*o). (2.30) 

Equation fl2.26p is convenient for several reasons. The solution of fl2.29p or (12.301) involves 
only a and is independent of the energy E. In other words, to obtain Mfi{E) one needs 
to solve these equations for a set of cr values irrespective to E. In fact only one of the 
equations is to be solved since l^^i) and |\E'2) are related to each other in a simple way. 
Furthermore, the dynamic calculation to be performed is independent of the final state 
channel /. 

Inserting and |\I^2) into (I2.26P one obtains Lfi{E,a). Then the Lorentz 

transform (12.260 is inverted, i.e. the integral equation ()2.24p is solved for Ffi[E, E'), 
and the singular integral in (I2.22p is calculated as 

POO 

/ {E-E' + t7]y'Ff,{E, E')dE' = -tnFf^iE, E) 

poo 

+ P (E - E')-^Ffi{E,E')dE' . (2.31) 

After having calculated the non-trivial term {(l)f{E)\{E — H + iri)~^0\'^Q) from (12.221) 
it is added in ( 12.2ip to the Born term to obtain the matrix element Mfi{E). [JJ 

2.2.1. Two-fragment break-up reactions. In case of two-body break-up reactions the 
calculation can also proceed in a different way. In fact this alternative method has been 
used in the various exclusive LIT applications discussed in section [HI 

The procedure is the same as described by (I2.18p - (l2.24p . but the Lorentz kernel is 
taken in the form ( 12.90 instead of ( 12.250 . Then one has 

Lfi{E,a) = {^a\^b), (2.32) 

with 

= {H-ay'd\^o), (2.33) 

\^,) = {H-a)-'\^f). (2.34) 

Obviously ^a,b are the solutions of the two inhomogeneous equations 

'H~aR-tai)\^a) = d\^o), (2.35) 



H-aR-taij = 10/). (2.36) 

In [21] the general and the two-fragment break-up methods have been used for 
calculating the d{e, e'p)n cross section. There were only small differences in the results 
(about 1%). 

II It is convenient to calculate the principal value integral entering here as a sum of ordinary integrals, 

PJZ ^dE' - ^'^E' - /l-A ^^^^dE' + /- , I^dE', where A is arbitrary, 

provided that Eth < E — A. 
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2.2.2. Case of fragments interacting via Coulomb potential. In order to take into 
account the long-range Coulomb interaction between the fragments one may proceed 
as follows. 

Let us denote with U'^ the average Coulomb potential between fragments, i.e. the 
Coulomb potential between the charges of the fragments concentrated in their centres 
of mass. Let be the continuum state given by the properly antisymmetrized 

product of the bound states of fragments times the 'plane-plus-ingoing wave' sub-state 
\ip~) of their relative motion in this average Coulomb potential. The latter is the wave 
function satisfying 

{frei + Uf -E)\ij-) = {), (2.37) 

where Trd is the kinetic energy of this relative motion. (.4.|<y9j ■*) may as well be the 
corresponding angular-momentum coupled state of this type). In analogy with fl2.20p . 
let us denote with ^) the properly antisymmetrized state 

|(^(-)) ^{H~ E)A\ipy^) ^ A{Vf - Uf)\v^f^) . (2.38) 

Then one may write (see equation (108) in chapter 5 of [19]) 

\^>-^{E)) = A\^y\E)) + {E-H- ^r^r'\^^f\E)) . (2.39) 

Because of this relationship one may use the same procedure as in sections 12.21 or 12.2.11 
with the replacements |0/) ^) and \(f)f) — * ^). However, one has to consider 

that, different from the case with no Coulomb interaction, the corresponding sub-states 
are no longer plane waves. This problem is easily overcome when the Coulomb 
interaction acts only between one pair of fragments since in this case is known 
exactly. 

2. 3. Comments on the application of the method to other kinds of processes 

With the integral transform approach one can also obtain the amplitudes of non- 
perturbative reactions of the general type, without calculating continuum states. In 
the non-perturbative case the procedure is rather similar to that described above for 
perturbation-induced transitions [1]. In fact in this case one has to calculate a matrix 
element like in fl2.18p . where, however, 0|\E'o) is replaced by \(t)i{E)), i.e. the initial state 
of the reaction, defined in a similar way as the state \(j)f{E)) in fl2.20p . The state |0j(i?)) 
normally corresponds to two colliding fragments. The key point is therefore again to 
calculate the form factor of fl2.23p type with the above mentioned replacement. The 
transformation with the Lorentz kernel may be employed to calculate this form factor. 
No such calculations have yet been performed. Since in this case the source term |0j(i?)) 
in the equations corresponding to fl2.27p and fl2.28p depends on i?, the corresponding 
dynamic calculation has to be done separately for each E value. 

Semi-inclusive processes, like e.g. (e, e'A^) reactions, are sometimes discussed in 
the framework of the so-called spectral function approximation [22]. Although the 
corresponding cross section is not an inclusive one, but only semi-inclusive, it may be 
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calculated with the procedure of the type presented in section I^TTl The calculation of the 
spectral function in the A-particle case is reduced to the corresponding calculation of a 
quantity of (12. ip type in the {A — l)-particle subsystem, with the state \Q) depending 
on the momentum of the knocked-out nucleon [23] . 

A non-perturbative amplitude of inelastic scattering of a fast projectile on a target 
in the framework of the Glauber approximation |23] is of the form f\Q), where \^ f) 
is the continuum final state of the target and \Q) includes its initial state and the 
profile function playing the role of the transition operator. This amplitude has the same 
structure as the perturbative amplitude of (12.181) and therefore can be calculated in a 
similar way. No such calculations are known. 



3. The LIT method: practical implementation 

3.1. Discrete and continuum spectra 

Separating in r(£'), defined in ( 12. ip . the discrete and continuum spectrum contributions, 
one may write 

^ {or - e„)2 + a\ {or - e)2 + o\ 

= LD{crR,aj) + LciaR,ai) (3.1) 

where we have introduced the excitation energy e = E — Eq. In particular e„ are 
the discrete excitation energies e„ = — Eq and r„ represent the corresponding 
contributions to the function r 

r„ = (Q|^„)(^„|Q')- (3.2) 
The second term in (13. ip is related to the continuum part of the spectrum with Cth 
representing the continuum threshold energy (e„ < Cth) and 

POO 

r(e)= / dE^{Q\^^){^^\Q')5iE^- E). (3.3) 

Jeth+Eo 

The aim of the LIT method is to obtain both r„ and r(e) solving (13. ip with 
L{(Jr, aj) as an input. It is expedient to use L calculated at a fixed o"/ in a range 
of aR values. One may also use cx/ = aj{aR). The first task is to get Lj:,{aR,aj) i.e. 
the discrete contributions to (13. ip . One way to get the contributions from the discrete 
levels is to calculate these levels explicitly to obtain and e„ and evaluate directly 
the overlaps (13. 2p and their contributions to the response function. Alternatively one 
can extract these contributions by calculating the transform L = (\E^|\E'') at a very small 
a I value for a range of aR values between zero and Cth- For fixed o"/, in the vicinity of a 
level e„, the transform L is dominated by the corresponding single term, 

LD{crR,cxi) c:^ " I 2 - (3-4) 

This aR dependence allows extracting e„ and r„. The procedure is convenient when one 
deals with few well separated levels. The accuracy in determining e„ and r„ may affect 
the continuous contribution r(e) to be obtained subsequently. 
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After calculating L{aR, aj) from fl2.13l) and subtracting the discrete contributions 
from it, the following integral equation of the first kind with the Lorentz kernel is to be 
solved to obtain r(e). 

f°° r(e) 

Lc{crR,o^i)= / de- — — 2- (3-5) 

To simplify the notation in the following we will omit the aj dependence and use L{(Tb) 
to indicate its continuum part only, i.e. Lc{cr^, unless specified explicitly. 

3.2. Role of the width of the Lorentzian kernel 

Let us comment on the choice of the o"/ value. The resulting response r(e) from an 
'exact' calculation of L and subsequent inversion should be independent of this choice. 
However, the rate of convergence in the calculation may depend on o"/. Let us show that 
the value of o"/ determines the relative error in r(e) manifesting in components with high 
frequencies. The LIT's of such high frequency components become very small because 
of the averaging out in the integral. Denoting with r,^exp(iz/e) the component of r(e) 
with frequency v and with its transform, one can show that if \i'\(Ji is less or about 
unity one has (o"//7r)|Lj^| ^ |r,^| while if |z^|o"/ ^ 1 one has 

— \Lu\ ^ ^ , ^ — ^TT^\^^- (3-6) 

When one solves (13.51) numerically one may adopt that the errors 5L in the transform 
and 5r in the response are related to each other in accordance with 



oo 



5r(e) 



5L{aR) = / de- — ^. (3.7) 

Let us denote 6Li, the component of 6L that corresponds to the component Jr,^ exp(ii/e) 
of 6r. Furthermore, let us use the notation 6Lq, Lq, 6rQ and tq for the components 
5Ly, L^, 5ry and with a low frequency. When \v\ai ^ 1, \Lu\ substantially 
decreases according to (13. 6p and it may occur that |(5L,^/(5Lo| ^ |L,^/Lo|. This gives 
3> |5ro/ro|. Therefore it is desirable that for all the frequencies z/ substantial to 
reproduce the main structure of r(e) the quantity |i^|c7/ would not be large. To this aim, 
in order to solve the integral equation (13. 5p a regularization is applied, which suppresses 
high frequencies (see section IXBl) . When the mentioned condition is fulfilled this ensures 
approximate stability of r(e) thus obtained and its closeness to the true solution. 

From the above point of view smaller u/ values are preferable. On the other hand, 
when a I decreases one approaches the scattering regime in (12.161) and (I2.17p . which 
makes it harder to obtain L with bound-state type methods. For example, suppose 
that those equations are solved in a subspace of localized functions of dimension 
(see section 13. 6p . If we denote = — Eq where E^ are the eigenvalues of the 
Hamiltonian matrix in this subspace then the approximate total transform (including 
both discrete and continuum contributions) has the structure 



N 
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Here the terms with lower values represent the discrete contributions in (13.11) while 
other terms correspond to a discretization of the continuum. When N increases the 
density of the latter type values increases as well providing a considerable number of 
them within a range (Jr — (Ji < ^ crR + o'j. When this occurs overlaps of peaks from 
separate contributions to (13.81) lead to a smooth approaching the true L as N tends 
to infinity. However, the smaller is cx/ the higher should be the value to represent 
the true L with a given accuracy. Therefore performing a calculation with too small a 
value of a J is not expedient. 

3.3. Role of the centroid of the Lorentzian kernel 

Let us comment on the range of values for which fl3.5p is to be solved. If one adopts 
that the spectrum r(e) to be obtained extends over the range Cf/^ < e < Cmax then it is 
reasonable to employ L^gr) approximately in the range eth — o-j < aR < Cmax + o"/- At 
such conditions complete information on r(e) contained in L{aR) is used. 

Suppose that r(e') does not change considerably in the range e — ctj < e' < e + aj 
and that the r values in this range are not much smaller than the values beyond this 
range. Then performing the integration in (13. 5p one may approximately take r out of 
the integral at the point e = aR. If e > Cth + crj this gives 

r(e) ~ —L{e). (3.9) 

71 

Thus the input transform itself looks similar to the output response. 
At large \aR\ the asymptotics 

LiaR) ^ ^ rie)de = ct^'{Q\Q'). (3.10) 

can be used to calculate the input transform (see the sum rule (12.51) ). However, one may 
prove that the asymptotics (13.101) is valid when and only when e^r(e) tends to zero as 
e tends to infinity. Moreover one can notice that at large \aR\ the regime (13.101) occurs 
when a^{Q\Q') is much higher then 7ir{aR)/aj while the regime (13.91) takes place in 
the opposite case. 

3.4. Remarks on particular cases 

Equation (12.41) corresponds to a non-relativistic calculation. It is implied that the center 
of mass subspace is separated out and states entering the response are defined in the 
relative motion subspace. The operators O and O' are 'internal' operators acting in the 
latter subspace. The corresponding matrix elements are taken between internal initial 
and final states with given momenta, with the momentum conservation delta function 
omitted. The energy E pertains to the relative motion subspace. In some cases, in order 
to take into account relativistic corrections one needs to consider responses defined 
in such a way, but with the energy conservation delta function that corresponds to 
relativistic kinematics. Up to some factors such responses may be expressed in the 
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same form as the non-relativistic responses fl2.4l) so that the present considerations are 
apphcable (see [25]). 

Performing closure to get f l2.8p we have assumed that the transition operators O 
and O' do not include dependence on the final-state energy E^. But in some cases, 
like e.g. in the electrodisintegration problem (see section 16.1.11) these operators are of 
the form GqOq,, where are form factors that depend on E^. In that case the 
response functions in the cross section may be written as 

^=5Z^"^-- + 5Z^-^/3(^"/3 + </3)' (3-11) 
a a<f3 

where Taa are functions of the form (12.41) with = 0' = 0^-, and r^p are again functions 
of that form with O = Oa, O' = Op. All the Tap are independent of E^ and are calculated 
as described above. Extra factors depending on E^ may also arise when the continuity 
equation is used to express some parts of the electric multipole operators of the current 
in terms of the density multipole operators [26]. In this case one may act in the same 
way. 

The operators O and O' may include a more complicated dependence on Er^. Then 
one may proceed as follows [27]. One replaces the operators 0{E^) and 0'{E^) with 
those 0{E) and 0'{E) depending on a fixed energy E rather than E^. The response r{E) 
is replaced with the arising subsidiary response r'(£^, and the latter is calculated as 
described above. If r{E) in some E range is of interest then r'{E,E) is to be found 
for a set of E values lying in the same range. After that r{E) can be obtained via 
interpolation as r'{E, E = E). 

3. 5. Inversion of the LIT 

As already pointed out, the main advantage of the LIT method to study reactions is 
that one avoids to solve the many-body scattering problem. One solves instead, with 
bound-state methods, equations of the form (12.161) . The knowledge of those solutions 
leads to the LIT of the function r of interest. A crucial part of the method is then 
the inversion of this integral transform. In section 13.11 we have seen how to obtain the 
discrete contributions to r, so one is left with the most challenging part represented 
by the solution of (13. 5p . The inversion of this integral transform has to be made 
with care, since it is unstable with respect to high frequency oscillations as explained 
in section 13. 2[ In this sense the LIT inversion problem belongs to the class of so-called 
ill-posed problems. 

In |5| the mathematical aspects of such problems are studied. It is also described 
how solutions of ill-posed problems can be obtained adopting a regularization scheme. 
Therefore, following [5] we have implemented a specific regularization in inverting the 
Lorentz transforms obtained in all the calculations reported in this review. This has 
led to very safe inversion results. Alternative inversion methods are discussed in |28j . 
They can be advantageous in case of response functions with more complex structures. 
Up to now, however, such complex structures have not been encountered in actual 
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LIT applications, since the various considered r(e) have normally (i) a rather simple 
structure, where essentially only a single peak of r(e) has to be resolved, or (ii) a more 
complicated structure, which however can be subdivided into a sum of simply structured 
responses, where the various LITs can be inverted separately. The second case (ii) has 
already been encountered in [3] and is also discussed in section I6.1.2[ 

The present 'standard' LIT inversion method consists in the following ansatz for 
the response function 

^ max 

r{e')= ^c„Xn(e',«i) (3.12) 

n=l 

with e' = e — Cf/i, where Cth is the threshold energy for the break-up into the continuum. 
The Xn are given functions with nonlinear parameters a^. A basis set frequently used 
for LIT inversions is 

X„(e,ai) = exp(- — ) . (3.13) 

n 

In addition also possible information on narrow levels could be incorporated easily into 
the set Xn- Substituting such an expansion into the right hand side of (13.51) (here too 
we omit to write the o"/ dependence) one obtains 

^ max 

H^r) = X] CnXn(o-i?,ai) , (3.14) 



n=l 



where 



oo 



Xn{crR,ai)=l de' ^"^^ ' "'^ . (3. 15) 

Jo (e' -(^r)^ + o-J 

For given the linear parameters c„ are determined from a least-square best fit of 

L{(Tfi) of equation (I3.14p to the calculated L{aR) of equation (I2.13P for a number of (Tr 

points much larger than Nmax- 

If one uses the basis set (13.131) . in performing the least-square fit of L{aR) one should 
vary the nonlinear parameter a2 over a rather large range. If no Coulomb interaction 
acts between fragments, with the proper ai the factor e"^ in (13.131) reproduces the energy 
dependence at threshold of a response function, or a form factor. This correct ai value 
can be obtained from the case when the short-range interaction between fragments is 
absent. Therefore, to have good quality results in the threshold region one may adopt 
this ai value. In the usual case when the lowest open channel is that of two charged 
fragments the factor in (13.131) may also be replaced with the corresponding Gamow 
factor to take into account the Coulomb interaction. Alternatively, one may include the 
parameter ai in the least-square fit. However, it should be noted that the least -square 
procedure may not be sensitive enough to a contribution from the threshold region to 
ensure the proper behaviour. 

For every value of Nmax the overall best fit is selected and then the procedure is 
repeated for N'^^^ = N^ax + 1 till a stability of the inverted response is obtained and 
taken as inversion result. A further increase of Nmax will eventually reach a point, where 
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the inversion becomes unstable leading typically to random oscillations. The reason is 
that L{aji) of equation (13.141) is not determined precisely enough so that a randomly 
oscillating r(e) leads to a better fit than the true response. If the accuracy in the 
determination of L{aji) from the dynamic equation is increased then one may include 
more basis functions in the expansion (13.141) . 

It is evident that the number of functions N^ax plays the role of a regularization 
parameter and has to be chosen within the above mentioned stability region. Normally 
such a stability region is reached without greater problems. However, in principle it 
can happen that one does not find a stable result. In such cases one can either try to 
improve the precision of the calculated L{aR) or use different basis sets, or enlarge the 
flexibility of the basis set (I3.13P by taking e.g. {x'^le, Oj), /5mm < /5 < Pmax}, 

Xi{e,ai) =e"'exp{-a2e), Xn(e, a*) = e^^ exp(-^), n>2 (3.16) 

with suitable values for Pmin and Pmax- 

In case that the inversion leads to a response that exhibits an unexpected structure 
it is useful to decrease the parameter aj in order to have a better resolution in the 
transform. 

We mention that according to what is said above any basis set may be used only 
up to a certain value of N^ax- Information which can be parametrized only via Xn with 
n > Njnax is lost anyway. Therefore it is important to work with a basis set, where a 
relatively small N^ax leads to a high-quality estimate of r(e). Varying the nonlinear 
parameters of Xn (and also /3 in case of Xn) equivalent to using many different basis 
sets. Usually it is no problem to find a proper set. 

The reliability of the inversion method outlined above has already been proved for 
various cases [31 EH EH]. A discussion about two additional interesting cases can be 
found in section [71 

3. 6. Calculation of the LIT via the eigenvalue method 

As we have seen in section [2] the LIT method can be used to reformulate a scattering 
problem as a Schrodinger-like equation with source terms which depend on the kind of 
reaction under consideration. These equations, which we will call the 'LIT equations' 
are essentially the same for any reaction, differing by the source term i.e. by their right 
hand side (see (I2.16p . (12.171) and (I2.36P ). In all cases the asymptotic boundary conditions 
are bound-state-like. Consequently the solutions of these equations can be found with 
similar methods as for the bound-state wave functions. Therefore the problem is much 
simpler than solving the Schrodinger equation for the continuum spectrum. Quantum 
mechanical bound-state problems may be solved in various ways. One can search for 
a direct solution of the dynamic equations (differential, integro-differential or integral) 
either in the coordinate or momentum representation. Alternatively, one can employ 
expansions over some basis set of localized functions. These expansion methods become 
more and more advantageous than the above mentioned ones as the number of particles 
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increases. Methods of all these types may be employed to calculate the localized LIT 
functions as well. In the following we discuss in greater detail the expansion methods 
to find the LIT functions. This is done both for the just mentioned reason and because 
they have been used for the vast majority of the LIT applications. (For a solution of 
the three-body LIT equation with the direct Faddeev-type approach see [30]). 

The truncation of the basis set converts the bound-state Schrodinger equation 
into a matrix eigenvalue problem and the LIT equations into a set of linear equations. 
These equations can be solved with various iteration methods and also Gauss type non- 
iteration ones. Such strategies have the drawback that one should solve these equations 
many times, as many as the number of ctr values that one needs for a proper inversion of 
the transform. In this section and in the following one we present two better strategies 
for calculating L{a). The first strategy, which we call the eigenvalue method and which 
we discuss here, involves the full diagonalization of the Hamiltonian matrix and expresses 
the LIT through its eigenvalues. This method is instructive from a theoretical point of 
view. 

Regardless of the reaction under consideration and the process that one wants 
to study the LIT method requires the calculation of the overlap (see f l2.13p - fl2.15p . 
or fl2:26|) - fl2:28|) . or fl232|) - fl2:34|) ) 

L{aR,ai) = (^1^') 

[H - Eq- aji + lai) [H - Eq- an - laj) 
where \Q') and \Q) contain the information about the kind of reaction one is considering. 
Seeking and as expansions over localized basis states, it is convenient to 
choose as basis states linear combinations of states that diagonalize the Hamiltonian 
matrix. We denote these combinations \(p^) and the eigenvalues . The index N is 
to remind that they both depend on A^. If the continuum starts at E = E^^ then at 
sufficiently high N the states \(Pn) having < Eth will represent approximately the 
bound states. The other states will gradually fill in the continuum as N increases. The 
expansions of our localized LIT functions read as 



ef: - Eq- an- taj 



Substituting fl3.18l) and fl3.19l) into fl3.17p yields the following expression for the LIT, 



and for \Q') = \Q) 



Li-n^-^) = \.^,._E^_,^y^,y (3-21) 
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From fl3.20p and fl3.2ip it is clear that L{aji, aj) is a sum of Lorentzians. The inversion 
of the LIT contributions from the states with < Eth gives the discrete part of a 
response function, whereas the inversion of the rest gives its continuum part. The 
spacing between the corresponding eigenvalues with > Et^ depends on and in 
a given energy region the density of these eigenvalues increases with N . (Since the 
extension of the basis states grows with A^, this resembles the increase of the density of 
states in a box, when its size increases.) For the reliability of the inversion procedure 
one needs to reach the regime when one has a sufficient number of levels within the 
(T/ extensions of the Lorentzians, as will be illustrated in section I7.3[ 

3. 7. Calculation of the LIT via the Lanczos algorithm 

In this section a second strategy to obtain the LIT is described. This utilizes the Lanczos 
algorithm [311 [32] to express the LIT as a continuous fraction. It turns out that for 
obtaining an accurate LIT only a relatively small number of Lanczos steps are needed. 
As the number of particles in the system under consideration increases the number of 
basis states grows up very rapidly and the Lanczos approach seems at present the only 
viable method to calculate the LIT. 

The motivation to use the Lanczos approach comes from the observation that from 
the computational point of view the calculation of the LIT is much more complicated 
and demanding than finding the ground state wave function of an A-particle system. 
In fact, to obtain the ground-state wave function one needs only to find the lowest 
eigenvector of the Hamiltonian matrix. On the contrary, as it is clear from (13.201) and 
(I3.2ip . the complete spectra of H over a wide energy range should be known to calculate 
the LIT. Therefore it is no surprise that the computational time and the memory needed 
to calculate L{a) have been the limiting factors in extending the LIT method for systems 
with more than four particles. It turns out that these obstacles can be overcome if the 
LIT method is reformulated using the Lanczos algorithm. Following reference [33], in 
this section it is shown how this can be done. 

To this end we assume that the source states \Q) and \Q') are real and rewrite the 
LIT in the following form 

L{a) = --Im ((g| ^\Q')\ . (3.22) 

A similar relation connects the response function r(e) to the Green's function 

rie) = --lm!.\imG{e + ir] + Eo)\ ■ G{z) = {Q\^—^\Q') , (3.23) 

provided that z = e + ir] is replaced by crji + iaj. This is not surprising since the 
properly normalized Lorentzian kernel is one of the representations of the 5-function 
and ai /tt L{(Tji) ?"(o"_r) for aj 0. In condensed matter calculations [Ml [35|, [36] 
the Lanczos algorithm has been applied to the calculation of the Green function with a 
small value of r], and its imaginary part has been interpreted as r(e) directly. This can 
be done if the spectrum is discrete (or discretized) and rj is sufficiently small. In our 
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case we have a genuine continuum problem and we want to avoid any discretization, 
therefore we calculate L{aR) in the same way, i.e. with finite o"/ using the Lanczos 
algorithm, but then we antitransform L{aR) in order to obtain r(e). 

This section is divided into two parts. First we discuss the application of the 
Lanczos method for a symmetric LIT, i.e. for the \Q') = \Q) case. This case is somewhat 
simpler then the general case discussed in the second part and can be used as a starting 
point for it. 

3.7.1. The Lanczos algorithm. The Lanczos algorithm [3T1 [32] is a technique to 
solve large, sparse, hermitian eigenvalue problems. The method involves partial 
tridiagonalizations of the matrix H without generating, however, intermediate full 
submatrices. Equally important, is that the extreme eigenvalues of H tend to emerge 
long before the tridiagonalization is complete. In order apply the Lanczos algorithm for 
calculating the ground state energy of an Hamiltonian matrix H one begins by choosing 
a starter vector |0o) (also called the Lanczos pivot), which must have a non-zero overlap 
with the ground state \ipo)- One then proceeds building the Lanczos orthonormal basis 
{\(f)i) ,i = 0, . . . ,n — 1} hj applying the Lanczos algorithm recursively 

= H\^,) - aM - , (3.24) 

where the Lanczos coefficients a,- and 6,- are defined as 



and 



\H\(Pi)- 6, =11 \w,) II; 6o = 0; , (3.25) 



,) = ^; |0_i) = O; (0.10,) = 5.,. (3.26) 



If one applies M Lanczos steps to an M x M matrix, then the original matrix is reduced 
into a tridiagonal form. As mentioned above, the power of the Lanczos algorithm lays 
in the fact that even after m « M steps the eigenvalues of the m x m tridiagonal 
Lanczos matrix are good approximations of the extreme eigenvalues of H. These 
eigenvalues converge rapidly with increasing number of Lanczos steps. Working with 
finite precision, it is necessary to recall that at each step there are round-off errors so 
that after a certain number of steps the Lanczos vectors loose their orthogonality. In 
order to obtain accurate results this loss has to remain small, otherwise one has to apply 
an additional reorthogonalization procedure. 

3.7.2. The symmetric case \Q') = \Q) . The application of the Lanczos algorithm to 
the symmetric LIT, 

L{a) = --lm\{Q\ ^|Q)| , (3.27) 

starts by choosing the normalized source vector 
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as the pivot for the Lanczos basis and setting z = Eq + +icrj. With the help of these 
definitions one can rewrite Lia) as 



L{a) = - — (QIQ)Im 



1 



z-H 

Therefore the LIT depends on the matrix element 

1 



(3.29) 



>o). (3.30) 

z — H 

This matrix element can be calculated applying Cramer's rule to the solution of the 
linear system (MJ [37] 

y^(^ - H)mnXnO = SmO , (3.31) 



(where Xno = (0n|^3^|0o)) which arises from the identity 

{z-H){z-Hr' = I, 
on the Lanczos basis i = 0, . . . ,n — 1}. Using Cramer's rule one gets 

det{Moo) 



(3.32) 
(3.33) 



M 



00 



(3.34) 



det{z - H) 

where 

/ 1 -h 

z — ai —62 
—62 -2 — 02 

V ; ; ; 

and the a„ and 6„ are the Lanczos coefficients. Defining Di as the matrix obtained by 
removing the first i rows and i columns from {z — H), one sees that Dq = {z — H), 

det{Moo) = det{Di) (3.35) 

and 

det{Do) = {z- ao)det{Di) - h\det{D2) . (3.36) 

Thus one has 

^ (3.37) 



2^00 



,2det{D2) 



z-ao- fer^gi(£,^) 

The recurrence relation f l3.36p . is valid for any submatrix Di. Substituting the 
appropriate expression for det{Di) one gets 

1 



^00 



z — ao 



(3.38) 



- "1 °-2det(D2) 



In this way one is able to write continued fraction containing the Lanczos 

coefficients and 6j, 

1 
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and thus also the LIT becomes a function of the Lanczos coefficients 

LW^-iij (3.40) 



^2~bl... 



This result illustrates how, using the Lanczos method, L{(7) is determined without 
solving the LIT equations, i.e. without inverting or diagonalizing the Hamiltonian 
matrix. 



3.7.3. The general case. Now we turn to the application of the Lanczos method to 
the general LIT case, i.e. when \Q) ^ \Q')- One can reformulate the problem in two 
different ways. 
Reformulation I 
The LIT 



L{a) = Im {Q\ 



z-H 



(3.41) 



is rewritten as a sum of four symmetric terms through the relation 



z-H 



ji{Q\ + {Q'\) ^ 

4 z — H 

-Jm+^{Q'\)^ 
4 z — H 

+ -m\-m\)^^ 

4 z — H 



^-^i\Q) + \Q'))-\m-{Q'\)^i 

4 z — H 



\Q) - \Q')) 



i\Q)+AQ')) 



(3.42) 



In this form one can apply the results of the previous section. 
Reformulation II 

Although quite elegant the previous reformulation requires four separate LIT 
calculations. In [33] a powerful alternative reformulation has been derived, which allows 
the calculation of all the non-diagonal LITs, just with one set of Lanczos coefficients. 
Writing the LIT in the form (see (12:251) ) 



Liz) 



1 



2iai 



{Q\ 



\Q'). 



(3.43) 



,z-H z* - H , 

choosing as previously the pivot (13.281) and assuming that after n steps the Lanczos 
vectors approximately form a complete basis, i.e. Yll=o 10*) (0*1 — o^i^ gets 

. ra-l 



L{a) 



2(Ti 



1=0 



Z-H 



\Q') 



z* - H 



\Q') 



2^{Q\(pi)lm{xio} 



(3.44) 



1=0 



where the pivot |0o) is defined as in (13.281) with \Q) — * \Q') and the matrix elements 

1 



z-H 



I0o) 



(3.45) 
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can be written as continued fractions of the Lanczos coefficients in a similar way as 
in (13.391) . One can notice that Imjxjo} = Imjxoj}, therefore in the following we describe 
a simple algorithm for a recursive calculation of all the xoi starting from xqo- 
Equation (13.311) can be generalized to 

^ ^ (^ H)rnn-^ni ^mi • (3.46) 



From this one can obtain all the XQi matrix elements for n 
element xqi is obtained solving the linear system 



/ z - ao 
-hi 


V : 



z — ai 
-62 





-62 
z - a2 



\ 



( a^oi \ 

3^21 

V '■■ J 



and m 



/0\ 

1 




The matrix 



(3.47) 



V ■■ J 



Using Cramer's rule one has 

det{Mm) 



det{z - H) 
where Mqi has the following form 

/ -61 
1 



(3.48) 



M, 



01 







- Oi 
-&2 



Z 





-62 
- 02 



Since 



one has 



V : 

det{Moi) 
det{z - H) 



(3.49) 



- hldet{D2) 



(3.50) 
(3.51) 

1 f3 52) 

{z - ao)det{Di) - h\det{D2) ' ' 

In section 13.7.21 it has been shown that for evaluating xqo, one has to rewrite it in terms 
of the ratio between det{D2) and det{Di). In an analogous way, Xqi must be written in 



bidet{D2) ; 
{z — aQ)det{Di 

bidet{D2) 



terms of the ratio 



detjPs) 
det{D2)- 



To this aim it is useful to rearrange the denominator of (13.521) 
in such a way that the ratio between det{D^) and det{D2) appears in the formula. This 
is achieved by using the following equation 



det{Di] 
In fact one has 



[z — ai 



(3.53) 



hi 

iz - a.Mz -aA-h^-iz- ar^hl'^^^ 
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Proceeding in an analogous way and using in addition the relation 
det{D2) ^{z- a2)det{D^) - hldet{Di) 

one obtains 



Xq2 = 



{z — ao)det{D-i) 


- bldet{D2) 


'{z - 02) 




ai){z - oo) 






bi 


1 f {z - ai){z - 






bi 




'{z - as) 




ai){z - ao) 


62 


^ b. 



-61)-— 



^^detiDs) 



bi 



bi 



1 [ { z - ai){z - ao) 
bi 



At this point it is clear that defining 



-bi 



bl 



z-a^- 



z-a4 -rr 



z - a„- 



"1^+3 



one has 



z - ao + 51(1) ' {z - ai)Aoo - 61 + Aoo5'(2) ' 

where Aqo — z — ao/bi is obtained from xoo- Now a new parameter 
{z - ai)Aoo - bl 



01 



can be obtained from xqi, leading to the evaluation of 

1 



{z - a2)Aoi - 62A00 + Aoi5'(3) ' 
which generates the parameter A02 

_ {z- a2)Aoi - 62A00 
^02 — 7 

and so on. Therefore one has 

1 

with 

Aoi 

where Aq-i = 1. 



z - ai)Xoi-i - biXoi-2 + Xoi-ig{i + 1) 
z — ai)Xoi-i — biXoi-2 



bi 



'i+l 



CONTENTS 



26 



4. Comments on related approaches 

4-1. Integral transform methods with other kernels 

A microscopic study of few-body responses with the integral transform approach has 
been suggested in [38111]. In [38] the Stieltjes kernel 

K{a,E) = -;^ (4.1) 

-C/ + (J 

has been employed to this purpose. The values of a are real here and —a does not 
belong to the spectrum. The corresponding transform 

^{a) = {m^\d\H + a)-'d\^,), (4.2) 

can be obtained with standard bound-state type methods. The inversion has not been 
considered. Instead it has been suggested to compare $(cr) with the corresponding 
transform of the experimental r{E). Such a procedure is inferior to the direct comparison 
between the experimental and theoretical r{E). In fact, when one passes from responses 
to transforms interesting physical effects, that can be visible only in certain limited 
ranges of E, can be obscured in the transform. In fact they can be spread out on a 
larger range by the integral with such a kernel. Moreover, experimental data in the 
whole range of E are required to perform the comparison and they may be missing. 
Sometimes the lack of data is compensated by theoretical results obtained within some 
model assumption. However, it may be difficult then to judge how appropriate these 
assumptions are, given the above mentioned spreading of that information by the 
integral. 

Besides the Stieltjes kernel, also the kernel 

'^>.^) = (£^. («) 

has been considered in In addition it has been shown how to calculate with this 
approach not only exclusive perturbation-induced reactions, but also general strong- 
interaction induced reactions. The necessity of an inversion of the transforms, in order 
to compare directly the responses with the experimental data, has also been underlined, 
since it is a necessary element for such applications. To this aim accurate inversions 
of the transforms with the kernels (14.11) . (14.31) have been performed in some model 
problems. However, it has been realized that a considerably higher accuracy in the 
input transforms than in the output responses was required. This is just the important 
difference with the case of the Lorentz kernel. 

In [6] , the realistic problem of the deuteron electrodisintegration has been addressed 
using the Stieltjes kernel (14. ip . The transforms (14. 2p have been calculated using 
standard routines for solving ordinary differential equations. It has been found that 
these transforms do not ensure sufficiently stable inversion. The conclusion is that the 
Lorentz transform is preferable even if the calculation of the Stieltjes transform may be 
somewhat easier. 
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The Laplace kernel 

K{a, E) = exp{-E/a) 



(4.4) 




(4.5) 



of r{E) defined in ( 12. ip . Monte Carlo techniques for calculating quantities of the form 
(14. 5 p have been studied starting from ^39j in connection with the determination of 
thermodynamic Green functionsliJ 

The problem of the inversion of the Monte Carlo Laplace transforms of temperature 
dependent condensed matter response functions has been addressed in |12], US]. The 
input transforms can be obtained with only a limited accuracy and the inversion is 
problematic. In |44J the transform (14. 5 p of a few-body (a-particle) response has been 
calculated using a known Green function Monte Carlo type technique. In [H] an attempt 
has been made to perform a comparison with experiment at the response level, inverting 
the Laplace transform of the response. However, information sufficient to judge on the 
quality of the inversion is not listed there. In subsequent papers the authors shifted to 
comparisons between theory and experiment at the level of the Laplace transform In [15] 
the Laplace transforms (14. 5p of electro disintegration responses of a number of nuclei have 
been compared with the transforms (12.61) with the kernel (14. 4p of experimental data. 

4-2. The method of moments 

The method of moments (MM) has been used to calculate responses in atomic physics 
(see e.g. j46l HZj). As an input one uses the moments 



They are calculated with the help of closure as explained in section [2711 The procedure is 
a recursive one. Comparing the MM with the LIT approach one finds a few similarities 
and one noticeable difference. The former consist in that also (14. 6 p can be considered 
as a kind of mapping depending on two parameters: in the LIT case an and cx/, in the 
MM case n and E. Moreover, in both cases one of them is kept constant {ai or E) and 
varies the other one {a^ or n). In this way one deals essentially with a mapping of E 
into cTfl or n. The important difference consists in the fact that the continuous variable 
E is mapped in our case to a parameter which can be varied continuously (o"/?) while in 
the MM case n can only be discrete. 

While positive results have been obtained with the MM approach, it has some 
limitations. The response is reconstructed from the moments using the classical 

% In [40 the authors refer to [41] inappropriately in relation to the present approach. Indeed, in [41] 
methods to calculate any many-body quantities like (|4.5p are not discussed. This paper is devoted 
exclusively to a procedure for reconstructing the Fourier or Laplace transform of the thermodynamic 
Green functions from some discrete Fourier coefficients. This procedure is based on the periodicity-in- 
time boundary condition of the thermodynamic Green functions at finite temperatures (the period is 
{kT)~^) and has no relevance to the present zero temperature problems. 




(4.6) 
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machinery. First of all one has to be sure that all used moments exist. Then the 
machinery is rather complicated and it includes solving nonlinear problems. Another 
shortcoming is related to the loss of accuracy in higher moments because of the necessity 
to calculate them recursively, while a rather high accuracy in its moments is required 
to reconstruct the detailed behaviour of r{E). 

Another interesting approach to calculate electron scattering response funnctions 
has been proposed in [H]. There a few 'cumulant expansions' have been suggested, that 
could restrict the number of moments necessary for reconstracting the energy spectrum. 
To our knowledge such an approach has not been investigated further. 



4-3. Response function from transform in the cr/ — ^ limit 
Using in (12. ip the Lorentzian representation of the delta function 

^(E,-E) = I lim,_ Im ^^^3^ = lim,_ ^ (^^ + (4-7) 

one realizes that r{E) may be calculated as follows 



r{E) = lim,,^or.,(^), r^^{E) = — (4.^ 

71 



where 



\^,^) = iH-E-taj)-'\Q), (4.9) 
\rj = {H-E-taj)-'\Q'), (4.10) 

Therefore r^jji^E) = {o-j/7T)L{E,ai). 

However, unfortunately the limit aj = cannot be realized in a simple way via 
an extrapolation of the o"/ dependent LIT. If, for example, the calculation of r„j[E) 
is performed solving (14.91) and (l4.1Up in a subspace of localized functions of a fixed 
dimension N for many o"/, then the extrapolation of its results to o"/ = would give a 
model response of the form 

N 

r{E) = Y,rnS{En-E). (4.11) 

n=l 

This limiting response in a given finite subspace is quite different from the smooth true 
response that corresponds to transitions to continuum. 

As already discussed in section 13.21 the smaller is aj the more difficult is to obtain 
I'^ai) and 1^^^^) and hence the response r„j{E), with the help of 'simple' bound-state 
type methods. Therefore we believe that this way of proceeding is more troublesome 
than the LIT method, which relies on the inversion of the transform. The latter gives an 
opportunity to perform a dynamic calculation with a higher aj value, which is certainly 
an easier task. 

In |l9] a procedure has been given to calculate the A = 3 response functions in 
the framework of the Faddeev description of dynamics. The equation of (14. 9 p form at 
(J/ — has been transformed to inhomogeneous integral equations of Faddeev type. 
Different from the case of aj finite, the solutions of such equations are not localized. No 
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inversion is required here. As a result the summation over contributions to a response 
from separate final states has been performed in a closed form, i.e. using the closure 
property, like in the LIT case. 

Here we should also mention calculations, made in the same spirit for A = 2 — 4, 
of the non-Born part of the T-matrix at complex energy with subsequent analytic 
continuation [50l [51] . 

4.4- Methods with Hamiltonian matrix diagonalization 

The response r{E) may also be obtained in the following way ||20j (compare 
also section 13. 6p . Consider a complete set of localized functions, and let \^^) and 
be eigenstates and eigenvalues of the Hamiltonian matrix in a subspace that spans 
the first states in the set. Then from (12. 7p one gets 

N 

^iv(a) = Y.^Q\v:!)K{a,e^.){v':\Q')- (4.12) 

V 

It is easy to show that if either (QIQ) or (Q'lQ') is finite and the kernel is uniformly 
bounded (A'(ct, < C) then $Ar — $(o") uniformly in cr as — oo. 

While the transform $7v tends to the true transform $ when A^ increases, the 
corresponding ri^{E) are quite different from each other at any A^. While the real 
r{E) is continuous and smooth, the corresponding to (I4.12p consists of transitions 
to quasi-levels (the same states as in (I3.18p - (I3.21I) ). 

N 

rN{E) = Y.^QWu){v':\Q')Ke^. - E). (4.13) 

V 

The transition amplitudes entering are rather chaotic, and as A^ — > cxd the tat do not 
tend to any limit in the normal sense. 

To get a reasonable approximation to the true r{E) from the approximate transform 
f l4.12p one should solve 

^n{o) = ^ K{a, E)r{E)dE, (4.14) 

imposing a regularization condition. This ensures a smooth output. In this way one 
comes to a continuous smooth response r|^r(-E) that does tend to r{E) as N ^ oo. When 
A^ increases the difference between the response r^(i?) and the response f l4.13p consists 
in components with higher and higher frequencies. 

In accordance with the discussion above it is expedient to apply this procedure 
with 'smoothing-type' kernels, like the Lorentz kernel (12. lip or the Gaussian kernel 
K{o-R, CT/, E) = exp[— (.E — o-rY/ct]] (for analogy with the Lorentz case the centroid and 
width of the Gaussian are indicated by cxji and aj, respectively). One may expect that 
$Ar((T) reaches stability at A^ values that provide sufficiently many levels in an energy 
range corresponding to the widths of the kernels. 

In some shell model studies (see e.g. [52]) the following procedure to calculate 
response functions has been adopted. The Hamiltonian has been diagonalized in a 
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large model space with the Lanczos algorithm and the model response f l4.13p has been 
obtained. Then this model response has been smoothed using the Lorentz or the 
Gaussian leading to the so-called Lanczos response (LCZR) smoothed function, and 
the result of fl4.12p form has been taken as the final response. However, such a quantity 
can reproduce the true response only up to the width of a smoothing function. On 
the other hand, the smaller is this width the larger dimension N of the model space is 
required to ensure convergence with respect to of such a smoothed response. 

In the light of our experience with the LIT method we recommend a better 
procedure: not to adopt the quantity fl4.12p as the final response r{E) but to extract 
the latter clS cl gularized solution to (14.141) . This would allow using the 'smoothing 
functions', or kernels, with larger widths and, correspondingly, employing shell model 
spaces of lower dimensions A^. Different from the above procedure, the resolution of the 
spectrum obtained in this way is better then the width of a smoothing function. This 
occurs because of the improvement in the results when one passes from a transform to 
a response. 

A numerical example which shows the comparison between the LCZR and that 
obtained by inversion of the LIT is given in section 17.21 



5. Solution of the relevant equations: the hyperspherical harmonics 
approach 

In section |2] it has been explained that the LIT method demands a good bound-state 
resolution technique that allows an accurate solution, not only for the ground state 
wave function of the target system, but also for the LIT functions (see (12.160 . (12.171) 
and ( I2.30p ). As was already stated in section 13.61 methods of expansion over basis 
sets turn out to be advantageous with increasing number of particles. In this section 
we shortly summarize one of them, namely the hyperspherical harmonics expansion 
method, since it has been used in most of the applications presented in this review (for 
LIT application to the No Core Shell Model (NCSM) approach see [SS]). 

We first introduce the hyperspherical formalism and define the HH basis on which 
the expansions are made. The two practical ways to accelerate the convergence rate of 



the HH expansion are described in the Appendix A and B. In particular, in [Appendix A 



we describe the correlated HH expansion method (CHH), and in Appendix B the method 
which makes use of an effective interaction defined in a subspace of the HH basis (EIHH) . 
Most of the results presented in this review have been obtained using one or the other 
of these two approaches. 

In the HH formalism, which has been successfully applied to study nuclear, atomic 
and molecular few-body systems [131 HH [13 HHl HTJ [IE], the Jacobi coordinates are 
replaced by a single length coordinate, the hyperradius, and a set of (3^4—4) hyperangles. 
The HH are the A-body generalization of the two-body spherical harmonics, and 
likewise depend only on the hyperangular (angular) coordinates in the hyperspherical 
(spherical) decomposition of the A-body (two-body) system. In general, the wave 
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function can be expanded in a series consisting of products of HH basis functions and 
hyperradial basis functions. This decomposition makes the HH expansion rather flexible 
in describing the A-body wave function and allows an accurate description of the wave 
functions at large distances. 



5.1. The hyperspherical coordinates 

To introduce the hyperspherical coordinates we start from the set composed by the 
centre of mass coordinate R = -j XliLi ^« normalized reversed order (A — 1) 

Jacobi coordinates (A represents the total number of particles) 



V2 



A-l 



A 



A-2 
A-l 



ri - ^ _ ^ (^2 + T-a H \-rA) 

r2 - ^ _ ^ (r-s + r4 H \- ta) 



riA-i = \ll{rA-i~rA) , (5.1) 



2 

which describe the positions of the j-th particle relative to the centre of mass of particles 
j + 1 to A. Each Jacobi coordinate ijj consists of a radial coordinate r]j and a pair of 
angular coordinates fjj = {9j, (pj). 

The (^4 — 1) Jacobi coordinates are then transformed into the hyperangular 
coordinates a2, ■ ■ ■ , ola-i through the relations 



sin q;„ = 77„/p„ , (5.2) 



where 



with Pq = 0. For n — A—1 one has 

p'^pU-\Y.^r,-r,f. (5.4) 

From this relation it is evident that the hyperradial coordinate p is symmetric with 
respect to permutations of the single particle coordinates. After the transformation 
described above the 3(^ — 1) internal coordinates for the A-particle system consist 
of one hyperradial coordinate p = Pa-i, and (3^ — 4) hyperangular coordinates 
^ = {^1) V2, • • • , VA-i, 0(2, as, ■ ■■ , aA-i}- 

By using hyperspherical coordinates one can write the Laplace operator for n Jacobi 
coordinates {n — 1 ... N), as a sum of two terms 

1 (9 (9 1- 

"--pl^-'dp/- dp^ pf-- ^'-'^ 
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The hyperspherical, or grand angular momentum operator K"^ of the n Jacobi 
coordinates can be expressed in terms of K^-i ^^"^ follows [5i 
,2 , 3n - 6 - (3n - 2) cos(2a„) d 



Kt 



+ 



dal 
1 



+ 



cos^ a, 



sin(2Q;„) 

• 2 n ' 

sm a„ 



(5.6) 



where we define Kf = i\. The angular momentum operator associated with these n 
coordinates is Ln = Ln-i + {Lq = 0). The operators K^, K^_^, L"^ and 
commute with each other. 



5.2. The hyperspherical harmonics basis 

The hyperspherical harmonics functions 3^[x„] are the eigenfunctions of the hyperangular 
operator (15. 6p . The explicit expression for the HH functions of the first n Jacobi 
coordinates is given by [55j 



{iimii2m2\L2M2){L2M2i3m3\L3M3) 

.mi,...,m„ 

n 

X {Ln-lMn M„) J]F,^,„^,(% 

n 



X 



J=2 



' 4.1 K I I 



sm an 



cosa,- 



\K 



j-i 



Pi'^"'^'''^-^'"^\cos(2a,)) 



(5.7) 



where Y^^ ^ are the spherical harmonics functions, pj"'*-* are the Jacobi polynomials 
and A/"^'* are normalization constants given by 



2(2/i + a + 6 + l)/i!r(/i + a + 6 + 1) 



(5.^ 



r(/i + a + i)r(/i + 6 + i) 

The symbol stands for the set of quantum numbers £i, -^^2, ■■■,Ln, /i2, •••,Ain 

and M„. The quantum numbers are given by 

Kj = 2/ij- + Kj_i + £j ; fxi=0 , (5.9) 

and the /Xj are non-negative integers. By construction, p^"3^[x„] is a harmonic 
polynomial of degree Kn- The HH function y[K„] is an eigenfunction of with 
eigenvalues Kn{Kn + 3n — 2). 

In view of equations (15. 3p and (15.41) it is evident that the harmonic oscillator (HO) 
Hamiltonian may be written in the form 



N 
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SAT - 1 9 

- + — + pM, (5.10) 



dp^ p dp p 

where fjj = rij/ro, p = p/tq where tq = (h/ {mto)Y^'^ is the oscillator length. It has the 
eigenvectors 

'^HO = RnMy[K] (5.11) 

with eigenvalues 

En = fuj ( ^^^~ ^\ r?j =!kJ (^ ^^^ ~ + + . (5.12) 

Therefore the HH i^T-quantum number can be associated with the quanta of excitations 
of the HO wave function, assuming no hyperradial excitation. 

When applying the HH method to the solution of the A-body problem one is facing 
two problems: 

(i) in general the HH basis states possess no special symmetry under particle 
permutations. For the three- and four-body problem, this is a minor problem since 
a direct (anti) symmetrization of the wave function is still feasible. For larger systems, 
direct symmetrization becomes impractical and there is a need for a more sophisticated 
symmetrization method [571 EH El] ; 

(ii) for strongly repulsive potentials the convergence rate of the HH basis turns out 
to be rather slow, therefore two methods have been suggested in order to 'accelerate' it. 
The first one (CHH) makes use of correlation functions, which correlate the pure HH 
basis in order to reflect the characteristics of the potential. In the second one (EIHH) 
the convergence is accelerated by an effective interaction, built from the bare potential 
by means of a similarity transformation. The description of these two methods can be 
found in Appendices A and B. 



6. Applications 

The LIT method has been applied to various perturbation induced reactions with nuclei. 
For the strong interacting nucleons the perturbation is represented by the electroweak 
probe (photons, electrons, neutrinos) and the cross sections contain the responses of 
the systems to the perturbation. In this section the formalisms for such reactions are 
outlined, putting in evidence how the cross sections are related to the response functions, 
which can be calculated with the LIT method. It is also discussed which are the relevant 
operators O and O' to be used for each reaction. 



6.1. Electromagnetic reactions 

The principal aim in studying of electromagnetic interactions with nuclei is to investigate 
the nuclear dynamics and to clarify the role of effective nuclear and subnuclear degrees 
of freedom. These are traditionally described by nucleons, mesons, and isobars and lead 
to two- and three-nucleon forces. In principle the effective degrees of freedom and the 



CONTENTS 



34 



corresponding forces should be related to the more fundamental quark-gluon dynamics 
of QCD. 

A great part of our present knowledge about the internal structure of nuclei has been 
obtained from electromagnetic reactions. The electromagnetic interaction has various 
advantages. It is a fundamental theory and it is weak enough, so that in most cases 
lowest-order perturbation theory is sufficient, allowing a simple interpretation of cross 
sections and other observables in terms of charge and current density matrix elements. 

In particular the study of the photoabsorption spectra has proved to be fundamental 
for the comprehension of the dynamics of a quantum system. Photoabsorption 
experiments on nuclear targets have been performed extensively since the dawn of 
nuclear physics until the beginning of the eighties. Later, however, this activity has 
slowed down somewhat. One of the possible reasons for this may have been the 
lack of theoretical progress, which made the comparison theory-experiment rather 
inconclusive. In relation to this, it is interesting to quote here a sentence from [60] 
that expresses this concern. The title of the paper is ' Photodisintegration of and 
^He\ In the conclusions the authors write: ^More theoretical work will be needed 
before the experimental results for the photodisintegration of three-body nuclei can test 
quantitatively the charge symmetry of the nuclear interaction. The present theoretical 
models are still rudimentary in that they do not include simultaneously a consistent 
treatment of the final state interactions...' and a list follows of other shortcomings in 
the calculations available at that time. The authors conclude that ^As a consequence 
of these problems the agreement between experimental and theoretical results is, in the 
mean, not notably good.' The progress achieved in recent years in the theory of few- 
body systems, and in particular the possibility opened by the LIT method of performing 
consistent treatments of initial and final states, also for systems with A > 3, has allowed 
a sort of revival of the field, with particular emphasis on the lightest nuclei, and has 
stimulated new experimental investigations with modern techniques. In particular, both 
inclusive and exclusive cross sections have been addressed. 

The information about nuclear dynamics that one can get from electron scattering 
reactions is complementary to that obtained by real photons and much richer. The 
reason is that, while the energy of the photon fixes the momentum transferred to the 
target, an electron can transfer energy and momentum independently. This allows to 
study the dynamics of the system at different wavelengths via the elastic or inelastic 
form factors (also called 'electron scattering response functions') that appear in the cross 
section. In about fifty years of experimental activity in electron scattering a wealth of 
data on targets all over the nuclear table has been collected at low as well as intermediate 
and high energies and momenta. Various laboratories are still active in this field, both 
in the United States (Jefferson Laboratory) and in Europe (Mainz, Bonn, Darmstadt). 
In Darmstadt a collider of electrons and beams of nuclei far from stability is planned. 
Considering both present and planned activities of these laboratories one can notice 
that they are concentrating on 'light' targets. A special role play the few-body systems 
and ^He, since they also serve as effective neutron targets. 
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6.1.1. Electron scattering reactions. In the one-photon-exchange approximation the 
inclusive electron scattering cross section can be written as [61] 



0"M 



%/?L(g,a;) + (|^ + tan2^)i?T(g,^) , (6.1) 



where gm is the Mott cross section, uj is the electron energy loss, q the magnitude of the 
momentum transferred by the electron, Q the electron scattering angle and = q^ — uj'^. 
The functions Rl and Rt are the electron scattering longitudinal and transverse response 
functions, respectively. The function R^ represents the response of the nucleus to the 
electron field via the nuclear charge and is given by 



) 



X 5{Ef -Eo + ^-u;), (6.2) 

where p{q, u) represents the charge density operator and Mt is the mass of the target 
nucleus. 

The results presented in section |8] are obtained with the following charge density 
operator 

A 

p(q,o;) = ^pf(q,.;)+p7(q,.;), (6.3) 



i=i 



where 



/3f(q,^)=e,e^'^-^ (6.4) 



= Gl{ql)^-±^ + Gliql)^-^ ^ I [GUi) + GliqDn,] , (6.6) 

A. = Gl,{ql)^-^ + GUql)^-^ ^ \ [GlMD + Gl{ql)T.,] . (6.7) 

In fl6.4l) - fl6.7p r, p, cr, and r are the nucleon position, momentum, spin and isospin 
operators, respectively, m is the nucleon mass, G^^/ are the proton/neutron Sachs 
form factors and G^\j are the isoscalar/isovector nucleon form factors. The two 
terms in fl6.5l) proportional to are the Darwin-Foldy (DF) and spin-orbit (SO) 
relativistic corrections to the main operator (16.41) (see e.g. [621 [63]). We refer to the 
main operator (16.41) as the non-relativistic one, although the dependence of the nucleon 
form factors on does not allow a non-relativistic interpretation. 

The operator p can also be conveniently rewritten in terms of the isoscalar and 
isovector charge nucleon form factors from (16.61) . 

p{q, u) = GUql)Ps{q) + GUqDpM)- (6-8) 
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As to the transverse response function, it represents the response of the nucleus to 
the electron field via the current density operator Jt transverse with respect to q 

Mq,^) = ^ df{^o\jUQ,uj)\^f){'^f\MQ,^)\'^o) 

X 6iEj -E, + ^-uj). (6.9) 

The current density in (16.91) consists first of all of the transverse components of a 
convection and a spin current 



with 



(6-11) 



' 2m ' 

3 

l-5^/i,^,xqe^''-^. (6.12) 

i 

The index (1) indicates that the operator is one-body. In order to satisfy the continuity 
equation, however, the momentum and/or isospin dependent n-body interaction V 
requires a n-body current density operator j^^^y This is called the interaction (or 
meson exchange) current. As is well known, the continuity equation is not sufficient to 
determine the ra-body current uniquely. Therefore in principle one needs a dynamical 
model for the nuclear potential which reveals the underlying interaction mechanism. 
Such a model can be supplied by the meson exchange picture of the nucleon-nucleon 
(NN) interaction. Even for a phenomenological potential a consistent two-body current 
can often be constructed [HH [65l [66] . 

If one wants to apply the LIT method to the calculation of the electron scattering 
response functions Rl,t one has to deal with an 'unpleasant' feature of the operators. 
In fact they depend on energy via the nucleon form factors. In section 13.41 it has been 
shown that the operators can be put in the form (13.111) and one can deal with the 
corresponding energy independent r^^ . 

6.1.2. Multipole Expansions. In order to apply the LIT method it is useful first to put 
in evidence the rotational quantum numbers J and M in each a/9-component of the 
response functions 

rf/5(g,^) = ^7^ E jtdf{^oiJo,Mo)\pUci)\^fiJ,M)) 

° Mo,J,M 

X (v[/^(J, M)|p^(q)|vI/o(Jo, Mo)) 5 (^Ef ~ Eo + ^ - u^j, (6.13) 
^api<l,^) = ^j^ E jtdf{MJo,Mo)\rj{^)\^f{J,M)) 

° Mo,J,M 

X (vI/^(J,M)g(q)|v[/o(Jo,Mo))5(^E/-Eo + ^-^) , (6.14) 
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where Jq is the total angular momentum of the ground state. Here / labels asymptotic 
quantum numbers other than J, M. The labels a, (3 refer to the isoscalar and isovector 
parts of the operators (see discussion about (13.111) ). Then one can expand the u- 
independent parts of the charge and current density operators (to simplify the notation 
we now omit the subscripts 'a' and '/3') as 

p(g)=47r^fp,„(g)y^:„(q), (6.15) 

jm 

jV(q)=47r (6-16) 

A=el,mag jm 

Here q = g~^q, and Yj"^ are electric and magnetic vector spherical harmonics |67j . 
The operators pjm and Tj^ are irreducible tensors of rank j. In (16.161) we set e = 
in the electric case and e = 1 in the magnetic case . With this choice the matrix 
elements of the Tj^ operators between the conventional angular momentum states are 



real. In Appendix C| it is shown that one can write 

2Jo + 



Att 



A-1T 

r^(?,^) = 77-T E E(2^+l)[^'^(^'^)]y' (6-18) 

A=el,mag Jj 

where the partial responses [r^yj and ir^y^ correspond to transitions induced by the 
multipoles pjrn and to states with a given J. There it is also shown that these 
partial responses do not depend on M and m. El 
The LIT of f leTT]) and fl618|l are 

^\l.on) = — --5^(2J+l)[L^(g,a^)]S, (6.19) 
^o'o + J- ^ 

L^{q,an) = —-- Yl Y^^J+lWiQ^^Ryf- (6.20) 

A=el,mag J] 

The partial transforms entering the right hand sides of (6.19) and (6.20) are calculated 
in subspaces with given J, M and parity from dynamic equations with sources depending 
on J and j. The calculation is M independent. To obtain the responses (16.171) 
and (16.181) from these transforms the integral equations of the form (13.11) are to be 
solved. Alternatively, one may invert the single terms separately and then sum up the 
corresponding partial contributions to the responses. 

+ The following notation is often found in the literature and will also be used in section 16.21 the 
longitudinal multipoles of order j are indicated with Cj, the transverse electric with Ej and the 
transverse magnetic with Mj. As an example: {r^Yj = rcj, {f'^Yj = '''Ej ii X = el and {r^Yj = 
if A = mag. Notice that in neutrino scattering (section 16. 2p also multipoles of longitudinal character 
appear in the response to the Zq field via the axial current density. This does not happen in electron 
scattering since charge conservation can be applied and such multipoles can be rewritten in terms of 
the Cj. 
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6.1.3. Photoah sorption reactions. The inclusive photoabsorption cross section is given 
by the following expression: 



where a is the fine structure constant, uj represents the photon energy and R{uj) is the 
response of the nucleus to the photon perturbation in unretarded dipole approximation, 
defined as 



Here |\I^o) and Eq are the nuclear ground state wave function and energy, \^ f) and Ef 
denote eigenstates and eigenvalues of the nuclear Hamiltonian if, and 



where A is the number of nucleons, rf is the third component of the isospin operator 
and Ti the coordinate of the ith particle in the centre of mass frame. 

As mentioned above, equations fl6.22p - fl6.23p represent the so-called unretarded 
dipole approximation to the total cross section. The approximations consist in the 
following: i) neglect of target recoil; ii) neglect of all electric multipoles Ej of order 
higher than j = 1; iii) after application of charge conservation, known as 'Siegert's 
theorem', reduction to the first order term in the Taylor expansion of the Bessel 
function of order j = 1 in Pjmil) (see [61]); iv) neglect of all magnetic multipoles 
Mj in the interaction between the photon and the nuclear currents. It is clear that 
approximations i)-iii) are justified for 'low'-energy photons. However, the limit of 
validity of approximation iv) needs to be investigated, because it is relevant in relation 
to the problem of nuclear currents and of their consistency with the potential, i.e. the 
problem, already mentioned in the introduction, of understanding the implicit degrees 
of freedom of nuclear dynamics. This can be understood considering the following: as 
already mentioned above in iii) the dipole operator ID is the long wavelength limit of 
the j = 1 term in Pjm{(l)- The longitudinal Ci operator (and not the transverse Ei 
multipole operator) enters the cross section, because one has used charge conservation 
connecting the current density to the charge density, thus one works with the so- 
called Siegert's operator ('Siegert theorem'). Therefore, while D is written in terms 
of single nucleon coordinates, information on many-body currents is included. In fact 
for energies where approximations ii) and iii) are valid, the knowledge of the form of 
the potential is sufficient to include all consistent currents in the dipole cross section, 
without knowing their formal expressions. However, Siegert's theorem does not apply 
to magnetic multipoles therefore if iv) is not valid one needs models for the many-body 
currents that can be tested in the comparison with experiments. 

Of course also the extension in energy of the validity of ii) and iii) has to be 
investigated. Some studies have been performed in the past on the two- and three- 
body systems [Ml EH], but the situation for heavier nuclei is still unclear. 




(6.21) 




(6.22) 




(6.23) 
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Regarding the exclusive reaction of a photon induced two-body break-up of the 
nucleus into a nucleon and a residual {A — 1) system, one can write the total exclusive 
cross section cr^^'^, within the approximations i)-iv) mentioned above, as 

a^^H = ^7r2ao;/c/i j dfi, | (vP" |£)|vI/o) P , (6.24) 

Mn,M(a-i) 

where n and k denote the reduced mass and the relative momentum of the fragments, 
respectively, Mjv and M(^a-i) the total angular momentum projections of the nucleon 
and the {A — 1) system, (a~i) final state continuum wave function of the minus 

type pertaining to the A^, {A — 1) channel [19] and E]y^(^A-i) is the energy of the final 
state 

En,{a-i) ^ ^ 2/i ~^ -^(^-1) ' (6.25) 

with -E'(A-i) representing the ground state energy of the fragment {A — 1). 

It is evident that R{lj) in (16.221) has the form (12. 4p . Therefore in order to calculate 
it, one has to invert the LIT fl21^ obtained by solving flZTB]) {O = 6' = D). In the 
exclusive case the matrix element in (16.241) can be calculated using the LIT method 
explained in section 12.21 again with O = D. 



6.2. Neutrino reactions 

The interest in inelastic neutrino reactions with nuclear targets stems from the role 
they play in major questions of contemporary physics. Such reactions are of central 
importance in various astrophysical phenomena, such as supernova explosions and the 
nucleosynthesis of the elements. 

Core collapse supernovae are widely accepted to be a neutrino driven explosion of a 
massive star. When the iron core of a massive star becomes gravitationally unstable it 
collapses until short-range nuclear forces halt the collapse and drive an outgoing shock 
through the outer layers of the core and the inner envelope. However, the shock loses 
energy through dissociation of iron nuclei and neutrino radiation, and gradually stalls, 
it becomes an accretion shock. It is believed, but to date not proved, that the shock 
is then revived as neutrinos emitted from the collapsed core (the proto-neutron star) 
deposit energy in the collapsing layers to overcome the infall and eventually reverse 
the flow to an outgoing shock which explodes the star. In order to revive the shock, 
the neutrinos must deposit about 1% of their energy in the matter behind the shock. 
The latter, which is assumed to be in thermodynamic equilibrium, is composed mainly 
of protons, neutrons, electrons, and light nuclei. In contrast to the fairly known cross 
sections of neutrinos with electrons and nucleons, the inelastic interaction of neutrinos 
with light nuclei such as ^He, ^He or cannot be accurately measured, and one must 
rely on theoretical calculations in order to evaluate the cross sections. 

The neutrinos migrating out of the proto-neutron star are in flavor equilibrium for 
most of their migration. The electron-neutrinos remain in equilibrium with matter for 
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a longer period than their heavy-flavor counterparts, due to the larger cross sections for 
scattering of electrons and because of charge current reactions. Thus the heavy-flavor 
neutrinos decouple from deeper within the star, where temperatures are higher. Typical 
calculations yield temperatures of ~ 10 MeV for /x- and r- neutrinos |69j, which is 
approximately twice the temperature of electron-neutrinos. Consequently, there is a 
considerable amount of z/^^,- with energies above 20 MeV that can dissociate the ^He, 
not to mention other nuclei, through neutral current reaction. 

The flux of neutrinos emitted in the collapse process is sufficiently large to initiate 
nucleosynthesis in the overlaying shells of heavy elements. Neutral reactions of alpha 
particle and neutrino in the inner helium shell are part of reaction sequences leading 
to the production of the rare A = 1 lithium and beryllium isotopes [70]. Thus, better 
understanding of the v-a reaction can lead to better prediction for the abundances of 
these elements. 

In the following we shall limit the discussion to neutral current neutrino-nuclei 
reactions. The extension to charged current reactions is straightforward, due to isospin 
symmetry. 

In the limit of small momentum transfer (compared to the Z-boson rest mass) 
and using the multipole expansion (see section [6.1.2[ including the footnote), the cross 
section can be written as [71] 

da 4G2 ... . . 

sm OdO X 




sm^^-|^cos^^) E + 



^^-cos2^5^2i?^.,,,.(g,c.) 



2g2 "~ 2 
J.j>o " ) 

where kfis the final momentum of the neutrino and 9 is the scattering angle. The rest of 
the notation is the same as for the electron scattering cross section (section 16.1.11) . The 
minus sign (— ) in (16.261) stands for neutrino scattering and the plus (+) for antineutrino. 
The functions RQQ,{q,uj) are the response functions with respect to the multipole 
transition operators O and O' (when O = O' we use the notation Rq = Rqq). At 
the supernova energy regime the neutrinos carry few tens of MeV. Therefore, it is both 
useful and instructive to apply the long wavelength approximation and restrict the sum 
in (I6.26P to the leading multipoles. In the limit g — > they are the Gamow-Teller Ef-^L^ 
and Fermi Cq operators which are g-independent and then the axial vector operators 
i?^, L^, M^, Lq and the vector CY , EY , L\ , which are linear in q . The main low-energy 
contribution usually comes from the Gamow-Teller and Fermi operators. However, for 
double closed shell nuclei such as the "^He nucleus these operators are suppressed and 
the leading contribution is linear in q. The complete current is composed of one-body 
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terms (the impulse approximation) and many-body terms. It turns out that in the long 
wavelength limit the leading axial two-body current contribute mainly to the Gamow- 
Teller term [72] . 

7. Numerical examples 

In this section we perform various numerical tests, both internal to the LIT method 
itself, and regarding the comparison with other approaches. We take as test 
observable the total photoabsorption cross section in unretarded dipole approximation 
(see section 16.1.31) of and ^He. In the first two sections two interesting test cases 
for the reliability of the inversion methods, as well as of the whole LIT approach, are 
presented. In section 17.21 we discuss the comparison between a response obtained via 
the LIT method and the corresponding LCZR of section WM . Finally, in section EH 
we present a comparison between the eigenvalue method and the use of the Lanczos 
algorithm to calculate the LIT (see sections 13.61 and 13. 7p . 

7.1. A test of the inversion algorithm on the deuteron 
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Figure 1. LIT of the deuteron photoabsorption cross section L^J of (|7.ip with cr/ = 10 
MeV. 

The deuteron total photoabsorption cross section is calculated in strict analogy 
to [6]. Following the notations there, one concentrates on the following LIT: 




Here only electric dipole transitions are taken into account {L = 1). Then one 
has to consider the following combinations for total angular momentum and angular 
momentum quantum numbers j and I of the outgoing neutron-proton {np) pair: (j, I) 
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= (0,1), (1,1), (2,1), (2,3). The inversion of leads to the unretarded dipole response 
function (16.221) and the corresponding cross section is given by (16.211) . The complex 
4>L=i,ji are obtained from the same radial differential equations as equations (27) and 
(28) of [6], but with a complex a = —(Jr + icx/ and taking for the excitation operator 9 
the dipole operator D of (16.231) . The radial differential equations can be solved in two 
different ways: i) direct numerical solution, as normally done for a two-nucleon problem, 
and ii) expansion of the solution in a complete set with a subsequent determination of 
the expansion coefficients. In this section we use method i), while method ii) is taken 
for the study in section 17. 2[ 

In figure [1] we show calculated with the AV14 NN potential [73] . One sees 
that the transform is dominated by a low-energy peak and that further structures 
are not present. An inversion of with the basis set (13.131) leads to a rather slow 
convergence in N^ax and thus we use the enlarged set (13.161) with ai = 3/2, Pmin = 1 
and Pmax = 2. In the past we did not observe a similar slow convergence for other LIT 
applications. This is caused by the steep low-energy rise of the cross section combined 
with a rather narrow peak at only about 2 MeV above threshold. In figure [2] we illustrate 
the convergence pattern of the cross section as a function of Nmax- One sees that the 
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Figure 2. Total deuteron photoabsorption cross section up to 50 MeV (a) and 15 
MeV (b) from inversion of LfJ with various -/Vmaa;-values (see (|3.12l) : 10 (short dashed), 
15 (dotted), 20 (long dashed), 25 (sohd), 26 (dash-dotted). 



peak height is steadily reduced with growing Nmax, while the cross section is quite stable 
at > 10 MeV. The reduction of the peak height ends with N^ax = 25; Nmax = 26 
leads to an almost identical result, and also for even higher Nmax the decrease of the 
peak height does not continue further. Therefore we consider the inversion result with 
Nmax = 25 as the final result for the cross section. 

As mentioned one can make a conventional calculation of the cross section by using 
the proper np scattering states. In figure [3] we show the result of such a calculation in 
comparison to that obtained via the LIT method, plotted in figure El One observes an 
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excellent agreement between the two calculations, showing the high precision that can 
be obtained with the LIT method. 




50 

CO [MeV] 



4 6 
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Figure 3. Total deuteron photoabsorption cross section up to 100 MeV (a) and 10 
MeV (b): LIT result (solid) and from calculation with explicit np final state wave 
function (x). 



7.2. Lanczos response vs. inverted LIT response 

In this section we discuss the comparison between a response obtained with the LIT 
method and the corresponding LCZR calculated as described in section 14.41 To this 
end we reconsider the deuteron total photoabsorption cross section (see section 17. ip . 
However, this time we use the non-local JISP6 potential [H]. This potential is 
parametrized in a harmonic oscillator (HO) basis and hence we also expand (pi^iji 
of (17.11) over the same HO basis, up to a maximal value of the HO quantum number 
Nho Ei'iid solve for with the Lanczos method, as explained in section 13.71 (Results 
for such a calculation with the JISP6 potential are already published in [75]). Here we 
discuss this case because it allows to demonstrate in a clear way the difference between 
a LIT inversion result and the LCZR of section U31 

In figure H] we show and the cross section with the JISP6 potential using 
expansions of 0^=1 with Nho = 150 and 2400. In both cases the maximal number of 
possible Lanczos steps, i.e. Nho, has been carried out. One observes that the increase 
of Nho has no visible effect on (relative differences are below 0.01% in the whole 
peak region). Thus it is not surprising that also the corresponding cross sections, are 
practically identical (inversions of L^J are made in a similar way as discussed for the 
AV14 case in section mi) . 

It is interesting to compare the obtained cross section with that of a LCZR 
approach. The results are illustrated in figure [5|, where for the LCZR a halfwidth aj of 
1 MeV is taken. With Nho=150 a very poor result is obtained, exhibiting very strong 
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Figure 4. Deuteron photoabsorption, (a): LIT of (|7.ip with ct/ = 10 MeV and 
(b): corresponding cross section (b). Upper panels: absolute results with Nho = 150 
(dashed) and Nho — 2400 (solid). Lower panels: relative difference of results with 
Nho = 150 and Nho = 2400 in %. 



unphysical oscillations. With A'^ijo=2400 the oscillations are considerably reduced, but 
they are still rather pronounced beyond 10 MeV. In addition the cross section is not 
reproduced correctly neither in the peak region nor at higher energies. One also notes 
a further deficiency: the cross section is non vanishing below the break-up threshold. 
If one decreases the halfwidth of the LCZR further to 0.5 and 0.25 MeV one does not 
gain much of an improvement, as shown in figure [5] (c),(d) for the A'^//o=2400 case. The 
oscillations become visible also at lower energies and their magnitude grows considerably. 
One would need to increase Nho much further in order to have a realistic result with a 
smooth and correct low-energy cross section. 

The oscillatory behaviour of the LCZR is due to the fact that only a tiny fraction 
of the Lanczos energy eigenvalues falls into the energy range of interest, while the 
vast majority of the about 10^ eigenvalues (there are four different (j, /) combinations 
(see text after (17.11) ) is located beyond 100 MeV. On the other hand we want to 
emphasize that Nho=^^^ is fully sufficient to have a converged LIT for aj of 10 MeV, 
which after inversion leads to a practically identical cross section as the A'^i7o=2400 case 
(see figure H]). This shows that via the LIT approach the information is processed in a 
much better way than with the LCZR. On the other hand for many-body calculations 
one finds cases in the literature (see e.g. [52|), where one can work also successfully 
with the LCZR. In these cases one has a huge number of states and Hamiltonians (e.g. 
effective s(i-shell interactions) that generate most of the Lanczos eigenvalues in the 
range of physical interest. 
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Figure 5. Deuteron total photoabsorption cross section, LIT result (thick solid) and 
various LCZR results: cr/ = 1 MeV with Nho = 150 (dash-dotted) and Nho = 2400 
(thin solid) in (a) and (b); Nho = 2400 and ct/ = 1 MeV (thin solid), cr/ = 0.5 MeV 
(dashed) and cr/ = 0.25 MeV (dotted) in (c) and (d). 
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Figure 6. ''He photoabsorption: LIT (ct/ 5 MeV) up to 100 MeV (a) and a zoom 
on the peak region (b) for ifo ~ 8, K^it — 9 (dashed) and for Kq = 14, Kljt = 15 
(sohd). 



7.3. A test of the LIT method on the a-particle 

In this section we illustrate a test of 'internal' consistency for the LIT method, when 
applied both to inclusive and exclusive cross sections. To this aim we calculate the 
total photoabsorption cross section of ^He in unretarded dipole approximation in two 
different ways, namely in the direct way as an inclusive reaction (section 12.11) leading 
to 0"*""^' and in an indirect way by summing the total cross sections of the '^He(7,p)^H 
and ^He(7,n)^He channels leading to a^^. The latter are determined by using the LIT 
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Figure 7. ^He photoabsorption, convergence of HH expansion for ct""^' (a) and a. 
(b): = 8, Klit = 9 (dashed) and Kq = 14, Klit = 15 (solid). 



formalism for exclusive reactions (see section [2.2.11 and 12.2.21) . Both approaches should 
yield identical results below the three-body break-up threshold, whereas beyond it 
one has cr*"^' > (j'^^ ., since also other reaction channels are open. Such a comparison 
between o"*"'^' and aT^^ has already been made [76j for a LIT calculation taking the MT- 
I/III potential [77j as NN interaction. A rather good agreement between a^"^' and 
has been obtained, but at very low energies non-negligible differences have been found. 
Since the low-energy region is particularly interesting, e.g. for reactions of astrophysical 
relevance, it is important to re-investigate this case. To this end we choose a very simple 
NN potential, namely the Volkov model [78], which is a soft-core central interaction. In 
addition we neglect the Coulomb force, thus one has = 2cr(^He(7,p)^H). The LIT 
calculation is carried out using HH expansions via the EIHH method (see Appendix 
B). In figure Owe show the LIT of the total cross section with ct/ = 5 MeV, where we 
take the following maximal grand-angular momentum quantum numbers K: i) = 8 
(ground state), Klit = 9 (LIT function) and ii) Kq = 14, Klit = 15. It is evident 
that the difference between the two LITs is very small (about 0.5 % in the peak region). 
Nonetheless one encounters important consequences for the inversion. For the case with 
Klit = 9 one obtains unstable inversions already for iVmax > 7, while the stabihty range 
extends to A/'max = 15 for the Klit = 15 case. Similar results are found for the various 
inversions which are necessary to determine a^^. 

In figure [7] we show results for the cross sections a*"''^' and o"!^'^ with the two different 
sets of K values. One observes that for a'^^'^\Ko = 8, Klit = 9) there is a steeper low- 
energy rise and a lower peak cross section than for a!^'^\KQ = 14, Klit = 15). On the 
contrary for a'^^ one finds a good agreement up to about 1 MeV, but somewhat different 
peak shapes. In figure M we compare al^^^ and a^^ for the same K value sets. It is 
evident that with the lower number of retained HH one does not have consistent results 




below the three-body break-up threshold of about 28 MeV, but with a higher number 
of them there is a much better picture. In order to judge the quality of the agreement 




CO [MeV] CO [MeV] 

Figure 9. Low-energy ''^He total photoabsorption, cr™'^' (a) and cr^^ (b) with Kq — 14, 
Kj^jT = 15. Inversion results with set (3.12) for aj — 5 MeV (solid), aj = 10 MeV 
(dotted) and from set (3.15) with ct/ = 5 MeV (dashed). In (b) also shown the Born 
result (dash-dotted). 

of the results we have also investigated the uncertainty in the inversions using set (13.131) 
with two different o"/ values (a/ = 5 and 10 MeV) and using the inversion set (13.161) 
with 0"/ = 5 MeV. For o"*"'^' with u < 28 MeV we find a rather energy independent error 
of about 0.08 mb. On the contrary the error on a^^ depends on u. It is close to zero 
below 23 MeV, while in the peak region it amounts up to about 0.08 mb. In figure [9] 
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we illustrate the situation below a; = 23 MeV. It is evident that one has very stable 
threshold results for cr^^, although the final state interaction (FSI) contribution is rather 
large as the comparison with the Born result shows. On the other hand cj^"'^' exhibits 
a rather large relative error. The origin of this rather different behaviour for cr^"'^' and 
cr^^ is due to the rather different LIT methods for inclusive and exclusive reactions. 
The main difference lies in the necessity to evaluate the principle value integral in (12.311) 
for the exclusive calculation. If the pole a.t E = E' is sufficiently far from the cross 
section peak, i.e. for the energies of figure [9], then small differences in the peak shape 
do not matter much in the integral. On the contrary, in case of the inclusive calculation 
one obtains one inversion result valid in the whole considered energy range and small 
differences in the peak shape propagate to the threshold region and can lead to rather 
large relative errors there. Of course in the present case one could also try to increase 
the precision of the calculation by further increasing Ko and Klit, which then could 
lead to similarly precise results as for the deuteron case of figure [31 

As mentioned above, the low-energy region is of particular interest. Therefore it is 
an important finding that a'^^ leads to very stable low-energy results: the convergence 
in K is rather good (see figure [7]) and the inversion error is quite small (see figure [9]). 

7.4- Eigenvalue versus Lanczos methods 

In the upper panel of figure [TO] we compare the inclusive LIT for the photodisintegration 
of "^He (see section I6.1.3P obtained with the eigenvalue method and with the Lanczos 
algorithm. The results have been obtained expanding the ^He ground state |\l/o) and 
the LIT functions using the HH basis and applying the EIHH method, described in 
Appendix B, to accelerate the convergence. For simplicity we have used the semirealistic 
central MT-I/III NN potential [77] . The ground state has been expanded on the set of 
four-body (anti) symmetrized hyperspherical harmonics (SYMHH) [5H] characterized by 
angular momentum L = 0, spin S = 0, isospin T = and = 0, up to grandangular 
momentum quantum number Kmax = 8. Since the dipole operator has j = 1 and no 
spin degrees of freedom, the LIT function has been expanded on the set of four-body 
SYMHH associated with the quantum numbers L = l, S = 0,T = 1 and = 0. The 
dipole operator only allows transition by one unit of K, therefore it has been sufficient 
to consider HH states with K < K^ax = 9. The LIT obtained with these values of 
Kmax have shown a sufficient convergence. The number of hyperradial basis functions 
has been 21 and the parameter aj has been set equal to 20 MeV. Since one difficulty in 
the Lanczos method is retaining the orthogonality of the basis vectors, we have checked 
this point and found that the overlap between the pivot and the lOOO'th vector is smaller 
then lO"'^. 

In figure [10 one sees that there is an excellent agreement between the results 
obtained in the two ways. In fact, up to 250 MeV, the maximal difference is lower 
than 5 ■ 10~^ (see lower panel). This picture deteriorates a bit when we decrease aj, but 
even for aj = 10 MeV the difference between the methods is better then 3 ■ 10~^ for the 
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Figure 10. Upper panel: ''He photoabsorption inclusive LITs calculated with the 
Hamiltonian inversion and with the Lanczos algorithm (150 Lanczos coefficients). The 
two curves are indistinguishable on this scale. In the lower panel the relative difference 
between them is plotted. 



same ctr range. Thus we can conclude that the two methods are essentially equivalent. 

In the following we discuss the main advantages of the Lanczos technique, which 
are i) its applicability to huge sparse matrices and ii) a considerable reduction in the 
CPU time. 

i) It seldom happens that one has to deal with a Hamiltonian matrix which is either 
sparse or written as an outer product. For such matrices the calculation of the product 
H\(t>) goes much faster then O(n^) (n is the dimension of the matrix), and there is a 
great incentive to use algorithms, like the Lanczos algorithm, that exploits this feature. 
In fact, the real power of the Lanczos formulation lies here: it opens up the way to 
work with a huge number of basis functions, which direct inversion methods or full 
diagonalization methods fail to address, ii) The reduction in CPU time results from 
the fast convergence of the lowest (and highest) eigenvalues, so that it is not necessary 
to perform all the Lanczos steps. In fact, checking the convergence of L{aji, aj) with 
respect to the number of Lanczos coefficients used in the continued fraction one finds 
the following: for (7^ — > oo only the first Lanczos step is necessary, because the entire 




Figure 11. Upper panel: convergence of the "'He photoabsorption inclusive LITs for 
a J — 10 MeV and for different numbers N of Lanczos coefficients. Lower panel: relative 
errors with respect to the converged Lanczos result (N=1500) 



continued fraction is dominated just by cr. Also in the low an region, as expected, the 
continued fraction converges rapidly, since for cr^ — the dominating contribution to 
the term l/{z — H) is given by the lowest eigenvalue of H, which is well approximated 
using only the first few Lanczos coefficients. It turns out that also in the intermediate 
an region a rather small number of Lanczos coefficients are sufficient. 

As an example in figures [TT] and [T2] we show the convergence patterns with respect 
to the number of Lanczos steps for the LIT of the inclusive photodisintegration of ^He 
described above. This time, the calculation has been carried out with all the SYMHH 
states up to K^ax = 15 and with 21 hyperradial states, i.e. a total of about 18, 000 basis 
functions. If written in matrix form the Hamiltonian would consume about 2.5 Gb of 
computer memory, and full diagonalization would be a long task. It can be seen that 
even for a sizable basis set the transforms converge very rapidly with only few hundreds 
Lanczos steps. As expected the convergence of the LIT is very good near threshold, 
then the number of Lanczos steps needed to achieve convergence grows with to a 
maximum value at about ~ 100 MeV and then, for higher values of convergence 
improves again. 
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To understand the effect of aj on the convergence pattern it is instructive to 
compare figures [TT] (a/ = 10 MeV) and [T2] (ct/ = 20 MeV). Roughly speaking, the 
convergence of the LIT improves by almost one order of magnitude when o"/ is increased 
from 10 to 20 MeV. This follows from the reduction in the sensitivity of the LIT to the 
exact position of the eigenvalues of H (see f l3.2ip ). when o"/ is large. In general for full 
tridiagonalization the Lanczos algorithm exhibits an 0(?t,^) time complexity. However, 
since the continued fractions converge already with a number of Lanczos coefficients 
much smaller than n, the overall time needed to get the Lorentz integral transform is 
proportional to ri^. In practice we have found that about 1000 Lanczos steps are enough, 
even for very large bases. 

In summary, reformulating the LIT method with the Lanczos algorithm represents 
an extremely powerful way to perform actual LIT calculations. All checks have shown 
that it is as accurate as any other method and it is much more efficient from the 
computational point of view. 
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8. Results 

In this section we present the resuhs of various LIT applications, to few-body systems 
of increasing number of nucleons. Those obtained in the two-body system mainly 
represent important tests of the LIT method. The others are results obtained within 
different potential models, ranging from semirealistic NN potentials (MT-I/III |77j . 
TN [9], AV4' and AV8' [79] (solid), MN [80]) to realistic ones (Paris [81], AV14 [73|, 
AV18 [H2|, BonnRA [83]). Also various versions of the 3N forces have been employed, 
namely the Urbana potentials UVIII [M] and UIX [85] and the Tucson-Melbourne 
(TM) [86]. 

8.1. Reactions with the two-body system 
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Figure 13. Longitudinal form factor RL{q,Enp) of tlie d{e,e'p)n reaction at 
q=440 MeV/c for Paris potential [81] of LIT calculation (dashed curves) and calculation 
with explicit np final state wave functions (full curves): contribution of CO and C2 
multipoles (a) and total results (b). 

We start the discussion reviewing the very first LIT application, i.e. the calculation 
of the deuteron longitudinal form factor Ri^q^uj) at g = 440 MeV/c. In [3J, where 
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Figure 14. Deuteron electrodisintegration (i(e, e'p)?T.. Real (left) and imaginary (right) 
parts of Coulomb transitions matrix elements of multipole order I (Born contribution 
subtracted) for a LIT calculation at g=440 MeV/c and Enp = 120 MeV with TN 
potential as a functions of the number of inversion basis functions N^ax'- diamonds 
{(J I — 5 MeV) , squares (ct/ ~ 20 MeV) ; also shown results for calculation with explicit 
np final state wave function (full curves) and deviations of ±1% from these results 
(dashed curves). 



the LIT has been originally proposed, it served as a test case for the applicability 
of the LIT approach. There, in order to have reliable inversion results it has been 
necessary to separate contributions from the Coulomb monopole (CO) and quadrupole 
(C2) transitions from the other multipole transitions (see also section 13751) . The former 
lead to a shoulder oi Rl close to threshold, while the latter exhibit a peak in the quasi- 
elastic region. In figure [I3] we show the results of this important LIT test (note: Enp is 
the kinetic energy of the outgoing np-pair in the centre of mass system). In the upper 
panel of the figure one sees that the CO and C2 multipoles exhibit a steep rise at the 
very threshold, leading to somewhat oscillating inversion results. In figure [T3](b) we 
present the comparison on the total Rl, where also the other Coulomb multipoles are 
taken into account. It is evident that there is a very good agreement between the LIT 
result and the result of a calculation with explicit np continuum state wave functions. 
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In [2T] a similar test study has been carried out for exclusive reactions. To this end 
the longitudinal part of the c?(e, e'N)N cross section has been investigated. In figure [T^ 
we show results for a few selected Coulomb transition matrix elements. In order to make 
the comparison more significant the Born contribution (see section 12.21) is subtracted. 
One sees that i) the calculation is essentially independent on the used crj-value and that 
ii) there is a very good convergence in Nmax in the inversion. Also shown are the results 
of a conventional calculation. Like for the inclusive test case of figure [T3] they agree very 
well with the LIT results. 

Among the two-body reactions considered within the LIT approach there exist 
also two interesting test cases concerning pion photoproduction. In [27j the KroU- 
Ruderman term has been considered and it has been shown that the results are the 
same as in a conventional calculation with explicit scattering states. A further study 
has been carried out in [89], where the proper implementation of A degrees freedom in 
pion photoproduction has been successfully investigated for a LIT calculation. 



8.2. Reactions with three-body systems 




Figure 15. Ratios i?L(with 3N force) /i?L(witliout 3N force) for triton at various 
momentum transfers; q values indicated in the figure in units of MeV/c. 

8.2.1. Electrons. For the three-nucleon systems first LIT calculations of a;) 
(g=250, 300 MeV/c) have been made using the Faddeev formalism and taking realistic 
NN potential models [30]. Like the deuteron cases of the previous section they have 
merely served as tests for the applicability of the LIT method. In fact an excellent 
agreement of the LIT results with those of a previous conventional Faddeev calculation 
has been obtained. Later a more systematic study of the effects of various theoretical 
ingredients on has been carried out with the LIT approach for a wider momentum 
transfer range using the CHH expansion technique [90l [25]. The dependence on the 
following theoretical ingredients have been investigated: NN potential models, 3N 
potential models, relativistic corrections to the non-relativistic charge operator (spin- 
orbit and Darwin-Foldy terms), reference frame dependence. Here we give a short 
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Figure 16. Comparison of theoretical and experimental Rl of ^He: AV18+UIX 
potentials (solid) and AV18 potential (dashed), experimental data from [87], up 
to three-body break-up threshold. Also shown theoretical results from [88] with 
AV18+UIX potentials (dotted). 

summary of the results of [90] . In general a rather small NN potential model dependence 
has been found, but in some cases there are also larger effects. These include the height 
of the quasi-elastic peak at lower q and the threshold behaviour. In figure[T5]we illustrate 
the effect of the three-nucleon force (3NF). It is typically between 5% and 10%, but 
reaches up to 15% for the low-energy responses, far beyond the quasi-elastic peak. The 
dependence on the 3NF model is very small. Relativistic corrections are unimportant 
at lower g, but cannot be neglected for q > 400 MeV/c. For momenta up to g = 500 
MeV/c the comparison of the theoretical results with experimental data is generally 
rather satisfying. The experimental data, however, are in most cases not precise enough 
to draw conclusions about the 3NF effect. A nice exception is the '^He low-energy 
response, where a 3NF proves to be necessary to obtain agreement with experiment (see 
upper panel of figure [T6!) . At higher q and low energy one finds a considerably higher 
Rl response in theory than in experiment (see lower panel of figure [16]). In figure [T6| 
for energies below the three-body break-up threshold, we also show results from a 
calculation with explicit three-nucleon final state wave functions [88], where the same 
theoretical ingredients as in our LIT calculation are taken. One sees that there are 
rather good agreements between the results of the two different calculations. 
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A study of the Rl reference frame dependence has been already initiated in (20], but 
in [25] has been worked out in more detail and with consideration of higher momentum 
transfers. To this end the i?^(a;fr, ^fr) of different reference frames has been calculated 
and then transformed in a relativistically correct way to the laboratory system according 
to 



RL{uj,q) 



^2 TTifr 

q 



(note the definition q) = R'1\lui,^, giab)) where and Mr are energy and mass 

of the target nucleus. In figure [17] we show results for four different frames [25], which are 




Figure 17. Frame dependence of Rl of ^He: LAB (dotted), B (dashed), AL (dash- 
dotted), ANB (sohd); experimental data from [HI] (squares), [S^] (triangles), [15] 
(circles). 



defined via the target nucleus initial momentum P^^: Py'^=0, = — q/2, P^^ = — q, 
pANB _ _3q^2. It is readily seen that one obtains more and more frame dependent 
results with growing q. As shown in [25] the frame dependence can be drastically reduced 
if one imposes quasi-elastic kinematics and takes the relativistic relative momentum 
between knocked-out nucleon and residual two-body system as input in the calculation. 
As seen in figure [181 this way one finds almost frame independent results, which show 
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a very good agreement with experimental data at q equal to 500 and 600 MeV/c, while at 
700 MeV/c the experimental peak height is considerably lower than the theoretical one. 
We would like to point out that the ANB result of figure fTTl lies within the band of almost 
frame independent results of figure [181 Thus one has to consider the ANB frame as the 
optimal frame for a non-relativistic calculation. In fact it leads to a proper description 
of the quasi-elastic peak region, with no need to assume quasi-elastic kinematics and 
thus its validity is not restricted to the peak region (further explanations why the ANB 
system is preferable are given in |25j). 




e),,(MeV) 



Figure 18. As in figure [iTl but considering two-body relativistic kinematics for the 
final state energy. 



8.2.2. Photons. Now we turn to the total three-nucleon photoabsorption cross 
sections. They have been studied with semirealistic NN potential models ^93] first, 
and subsequently also with realistic nuclear forces (realistic NN potential + 3NF) |95j . 
In both cases CHH expansions have been used. In figure [11] we report the results of the 
first calculation of the total photoabsorption cross section with realistic forces [95]. The 
cross section is rather independent from a complete nuclear force model with NN and 3N 
forces. The role of 3NF is not negligible, they lead to a reduction of the low-energy peak 




Figure 19. Total photoabsorption cross section of (upper panel) and '^He 
(lower panel) with AV14 (dotted), AV14+UVIII (solid) and BonnRA+TM (dashed) 
potentials; the shaded area represents the data from [6^ while the black dots are the 
data from [53] 



by about 7% and to similarly strong increases at higher energies. One finds a rather 
good agreement of the theoretical results with experimental data. The experimental 
error bars, however, are quite large. In [29] the AV18+UIX LIT results, which are very 
similar to the AV14+UVIII and BonnRA+TM results of figure fT9l have been successfully 
benchmarked with the Faddeev results of the Bochum-Cracow collaboration. In [29] also 
a check of the validity of the unretarded dipole approximation has been made. It has 
been found that retardation and higher multipole contributions are quite small up to 
moderate photon energies (below 2% for u; < 50 MeV). 

A high-precision comparison of CHH- and EIHH-LIT results with the AV18+UIX 
nuclear force model has been made in |96] . Finally we would like to mention that various 
photonuclear sum rules (polarizability, bremsstrahlungs, and TRK sum rules) have been 
studied in pll[95]. 

8.2.3. Neutrinos. Next we consider weak reactions with three-body nuclei. The 
calculation of the cross section for neutral current neutrino scattering of '^H and ^He 
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Table 1. Averaged neutrino- and antineutrino-'^H and ■^He neutral-current inclusive 
inelastic cross sections per nucleon {A — 3) (cr)T„ — jj {<^i^ + '^t7)t„ (left columns), 
and energy transfer cross sections, {u)a)T^ = {ujct^ + uj(t-)t„ (right columns), as 
a function of the neutrino temperature T^, in units of 10^*^ cm^ and 10^**^ MeVcm^ 
respectively. 



in the supernova scenario with reahstic nuclear Hamiltonian and MEC currents has 
been presented in [97]. In the supernova scenario it is customary to assume that the 
neutrinos are in thermal equilibrium. Consequently, their spectra can be approximated 
by the Fermi-Dirac distribution with a characteristic temperature T^, and with zero 
chemical potential 

/(^^'^) = |3^^7ItI' ^'-'^ 

where k is the neutrino momentum and A^~^ = 2 ^^-^(— l)""'"^/n'^ is a normalization 
factor. In order to estimate the mean neutrino interaction with the nuclear soap 
the calculated cross sections must be averaged over incident energy and angle. The 
quantities of interest are the temperature-averaged cross sections and energy transfer 
cross sections: 

poo p 1 

Wt. = dcu dhf{n,ki) — , (8.3) 

{uj(t)t,= I duj dki f{T^,ki) u — , (8.4) 

where fcjj are the initial and final neutrino energy, uj = ki — kf is the energy transfer, 
and Uth denotes the threshold energy of the break-up reaction. In Table [H we present 
the results of [HZj for the averaged neutrino and anti-neutrino neutral-current inclusive 
inelastic cross sections and energy transfer cross sections as a function of the neutrino 
temperature. The difference between the '^H and the ^He cross sections reflects the 
difference in thresholds between the two nuclei. The mirror symmetry between both 
nuclei is restored with higher neutrino energy. The leading contributions to the cross 
section are the axial E^, Mf and E2 multipoles. The relative importance of these 
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multipoles varies as a function of the momentum transfer, and thus as a function of 
the neutrino temperature. In comparison to inelastic excitations of ^He studied in |98j . 
the cross sections (the mean values of and ^He) are about a factor 20 and 10 times 
larger at temperatures of 4 MeV and 6 MeV, respectively and the energy transfer cross 
sections are 8 and 2 times larger. 

At low momentum transfer, the Gamow-Teller operator dominates for the cross 
section, consequently the MEC have a large effect of about 16% at a temperature 
of 1 MeV. At higher momentum transfer, higher-order multipoles start to play an 
important role. Due to spatial symmetry, the MEC contribution to these multipoles is 
small and the overall effect of MEC decreases rapidly to less than 2% for temperatures 
above 4 MeV. While not directly important here, the asymmetry between the scattering 
of neutrinos and anti-neutrinos increases with temperature: the difference in the energy 
transfer cross sections grows gradually from 3% for a neutrino temperature of 3 MeV 
to more than 50% for a temperature of 10 MeV. Finally, the cutoff dependence of these 
observables is less than 2% for 1 MeV and less than 1% for higher temperatures. The 
estimated precision of these predicted cross sections is a few percent, which also includes 
estimates of the numerical accuracy. 

8. 3. Reactions with the four-body system 

8.3.1. Electrons. The longitudinal form factor RL{q,uj) of ^He has been calculated 
in |101] at 300 MeV/c < q < 500 MeV/c with the semirealistic TN potential via the 
LIT approach and using the CHH expansion technique. Already before, calculations 
with a realistic nuclear force had been carried out using the Green function Monte 
Carlo approach for the Laplace transform of the response [H]. However, the inversion 
of this integral transform is problematic (see discussion in section l^?T|) . In [101] a rather 
good agreement of the LIT results with experimental data has been found, as can be 
seen in figure [201 The figure also contains LIT results from [23], where Rl has been 
determined calculating the spectral function S{E, k) (see section 12. 3p of the residual 
three-nucleon system in the "^He(e, e'p)X reaction. In the calculation with a spectral 
function the FSI of the outgoing proton with the residual three-nucleon system is not 
taken into account. The figure illustrates that the proton FSI is quite important. It 
leads to shifts of the quasi-elastic peak positions by about 15 MeV towards lower energy 
and reduces the peak heights quite considerably (30%, 20% and 10% at g = 300, 400, 
and 500 MeV/c respectively). Rather strong FSI effects are also found close to the 
break-up threshold. 

Very recently also Rriq, t^) has been studied [26] with the MT-I/III potential, using 
EIHH expansions. Besides the non-relativistic one-body currents (lA), a consistent 
meson exchange current (MEC) is taken into account explicitly, (with an approximation 
which, however, is shown to be very good at lower momentum transfer). In figure [21] 
we illustrate Rt for different momentum transfers in comparison to the available 
experimental data. The spin current strongly dominates the response function at higher 
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Figure 20. Rl of ''He with the TN potential: plane wave impulse approximation result 
via use of spectral function (dashed curves) and full results (solid curves); experimental 
data from [99] (diamonds) and |100j (circles). 

momentum transfer (here not shown for q = 300 MeV/c, for more details see t26j), while 
at g = 100 MeV/c the convection current gives a large contribution. One also notes 
quite an important MEG contribution: at energies of about 100 MeV MEG lead to an 
enhancement by a factor of four with respect to the lA. With increasing q the enhancing 
effect of the MEG with respect to the lA is reduced dramatically. Unfortunately, there 
are no experimental data for the lower two g-values. At g = 300 MeV/c one finds a fair 
agreement with the data at lower and higher energies, but the quasi-elastic peak height 
is somewhat underestimated. 

For the four-body system two electron induced exclusive reactions have been 
considered with the LIT method (GHH expansion technique) and using semirealistic 
NN potentials: the longitudinal parts of the cross sections of ^He(e, e'p)^H, close to 
quasi-elastic kinematics [103] . and of ^He(e, e'(i)^H, at rather low u and 250 MeV/c 
< q < 450 MeV/c |104j . For ^He(e,e'p)^H results on the FSI effects and a comparison 
with experimental data are shown in figures [22] and [231 In figure [22] the considered 
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Figure 21. Rt of ''He with tiie MT-I/III potential: spin current only (dotted), 
lA with spin and convection currents (dashed), full calculation with the inclusion 
of a consistent two-body current (solid); experimental data from [SS] (open circles) 
and [100] (solid circles). 



kinematics, labelled with numbers from 1 to 9, are illustrated in the q-uj plane. 
In figure [23] the effects of the FSI for the nine cases are shown. One sees that FSI 
reduces/increases the cross section below/above q = 500 MeV/c leading to an improved 
agreement with experiment only for q below 500 MeV/c. Presumably for higher q 
an inclusion of relativistic effects might be important and possibly lead to a better 
agreement with experiment. 

For the (e, e'd) reaction of figure [24] one finds an enormous reduction of the cross 
section due to FSI, which, however, is not sufficient to describe the rather small 
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Figure 22. *He(e,e'p)^H reaction: position of the various kinematics of figure [23] 
with respect to the u = g^/(2mp) ridge. Shaded area represents the region below the 
break-up threshold. 
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Figure 23. ^He(e, e'p)^H reaction: Percentage deviation of theoretical differential 
cross section from the experimental values of |102| : without FSI (open circles) with 
inclusion of full FSI (solid circles). 



experimental cross section at lower q. As discussed in |104j a consideration of the tensor 
force could be particularly important for the FSI and possibly lead to an improved 
agreement with experiment. 

8.3.2. Photons. The ^He photodisintegration has a rather long history and many 
different experiments have been carried out over the years (see e.g. [TS]). Unfortunately 
the results of the various experiments on the total photoabsorption cross section 
differ quite substantially in the size of the giant resonance peak. About ten years 
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Figure 24. '*He(e, e'd)d reaction, differential cross section at a kinetic energy of the 
outgoing dd-pair of 35 MeV and averaged over a missing momentum range between 
100 and 150 MeV as a function of (four-momentum transfer squared): without FSI 
(dashed) , with FSI (solid) . Experimental data from [105] . 

ago theory and experiment [1091 1110] seemed to converge to a rather low peak 

cross section. However, in the LIT calculation of [113j a much more pronounced 
giant resonance peak was found. The calculation had been carried out in the 
unretarded dipole approximation with semirealistic NN potentials models (TN and 
MT-I/III potentials) using CHH expansions. We would like to point out that the 
unretarded dipole approximation should work even for higher energies than for the 
total photoabsorption cross sections of two- and three-body nuclei, since the ^He root 
mean square radius is smaller than the radii of ^H, and ^He. An improved LIT 
calculation has been later pursued using EIHH expansions [114j . Somewhat modified 
results have been found, but the pronounced giant dipole peak cross section has been 
confirmed. The relative importance of the considered reaction and the contradiction of 
theoretical and experimental results have led to a renewed interest in the low-energy 
^He photodisintegration. New experimental data |108l I112j as well as LIT calculations 
for the exclusive ^He(7,p)^H and ^He(7, n)^He reactions (CHH expansions) [76] and a 
first calculation for the total ^He photoabsorption cross section with a realistic nuclear 
force (EIHH expansions) |119j has become available thanks to the LIT method. The 
present situation for the total cross section is summarized in figure [251 One observes a 
rather good agreement of the theoretical calculation with the data of [108j . while there 
is a large discrepancy with the data of [112] . One further notes that the effect of the 
three-nucleon force is quite large, it reduces the giant dipole peak height by about 10% 
and increases the high-energy cross section quite considerably (35% at pion threshold). 
The comparison of the full theoretical result with the data of [106] is quite satisfying, 
although the data shows an even higher peak cross section. 

Very recently a LIT calculation using chiral nuclear forces (NN+SN) [120] has been 
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Figure 25. Total "^Hc photoabsorption cross section: (a) MT-I/III (dotted), AV18 
(dashed), AV18+UIX (solid). Experimental data from |106| . (b) as (a) but with 
upper/lower bounds included (estimate of theoretical error of EIHH convergence); 
experimental cross sections: dotted box [107! , squares |108j . open circles |112| . The 
area between dotted hnes has been obtained combining data from |109j and [llOj . 



carried out. The total ^He cross section results are very similar to the AV18+UIX 
results of figure [251 Calculations with non-local nuclear force models, where a 3NF 
is eliminated by construction (JISP6 and UCOM [121] NN potentials) have also been 
performed with the LIT method [TSl 1122] . They confirm the pronounced giant dipole 
cross section, but with peak positions at somewhat higher energies. 

In |123] the integrated ^He total photoabsorption strength has been considered. In 
particular the polarization, bremsstrahlung and TRK sum rules with the AV18+UIX 
nuclear force have been investigated. With their help proton-proton and neutron- 
neutron distances have been determined and a small deviation from the tetrahedral 
simmetry of the spatial configuration of ^He has been found. 

The exclusive reactions ^He(7, A^)^X have been calculated with the semirealistic 
MT-I/III potentials. The sum of the two channel cross sections can be compared with 
the total ^He photoabsorption cross section. Below the three-body break-up threshold 
both results should be identical, since other channels are not open (^He(7, d)d cross 
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Figure 26. More-body break-up cross section in ''He photodisintegration: total 
theoretical result of three- and four-body break-up (dotted); experimental results 
for *He(7, pn)(i (full symbols) and "'He(7,2p2n) (open symbols), [115] (upward 
triangles), |116j (squares), |117j (circles), (llSj (downward triangles). 

section is zero in unretarded dipole approximation). On the other hand at higher 
energies the difference between the total ^He(7) and the ^He(7,p)^H plus ^He(7,n)^He 
cross sections leads to an estimate for the more-body break-up cross section. The 
following results have been found in [76] . Below the three-body break-up threshold the 
agreement between the cross sections of the two kinds of calculations is rather good, but 
not yet completely satisfying. Therefore in section 17.31 we have reconsidered a similar 
consistency check of the LIT calculation. There it is shown that the HH expansion has 
to be carried out to a somewhat higher value of the grand angular HH quantum number 
K than has been done in [76] in order to find a satisfying agreement between both 
calculations. Concerning the more-body break-up cross section, in figure [26] we show 
the sum of three- and four-body channel cross sections obtained as described above. 
One sees that there is a fair agreement between experimental and theoretical results. 

8.3.3. Neutrinos. The inelastic neutral current neutrino scattering from ^He has been 
calculated in [124j in the impulse approximation using the AV8' NN potential. Later 
in [98] the full nuclear Hamiltonian with the AV18+UIX potential has been used and 
MEG have been included. 

In the supernova scenario one has to consider neutrinos with up to about 60 MeV. 
Usually, the leading contributions in weak nuclear processes are the Gamow-Teller and 
the Fermi operators. Due to the total angular momentum and spin structure of the 
^He nucleus, they are both strongly suppressed. In fact, the Gamow-Teller operator 
contributes only due to the small P and D-wave components of the ground state 
wave function. In addition, '^He is an almost pure zero-isospin state [1251 1126] , hence 
the Fermi operator vanishes. Therefore, the leading contributions to the inelastic cross 
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Table 2. Temperature averaged neutral current inclusive inelastic cross section per 
nucleon (in 10^^^ cm^) as a function of neutrino temperature (in MeV). 

section are due to the axial vector operators E^-, Mj^, Lq , and the vector operators 
CX ,E\ ,L\ (the latter are all proportional to each other due to the Siegert theorem). 
For the neutrino energies considered here it is sufficient to retain contributions up to 
0{q^) in the multipole expansion |124j . In figure [271 we present the convergence in Kmax 
for these multipoles. It can be seen that the EIHH method results in a rapid convergence 
of the LIT calculation to a sub-percentage accuracy level. 




Figure 27. LIT convergence in Kmax for the leading multipoles. 

In Table [2] we present the temperature averaged total neutral current inelastic 
cross section as a function of the neutrino temperature for the AV8', AV18, and the 
AV18+UIX nuclear Hamiltonians and for the AV18+UIX Hamiltonian adding a MEC 
derived from a chiral EFT Lagrangian [127] . From the table it can be seen that the low- 
energy cross section is rather sensitive to details of the nuclear force model (the effect of 
3NF is about 30%). This sensitivity gradually decreases with growing energy. In contrast 
the effect of MEC is rather small. Similar behaviour has also been observed for the hep 
process [128j. The small contribution of the MEC can be understood in the following 
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way. The symmetry of the vertices in this low-energy approximation dictate a symmetry 
between the two nucleons interacting via the meson exchange. The leading one-body 
multipoles have negative parity, as a result the MEC contributions to them is small. In 
comparison the MEC correction to the Gamow-Teller is of the same magnitude as the 
one-body current, however both terms become marginal with increasing momentum 
transfer. Although presented for the neutral current, these arguments hold true also 
for the charged currents since the response functions are related by isospin rotation. 
These results are the first fully microscopic study of z/ — a reactions, using a state of the 
art nuclear Hamiltonian and including MEC. The overall accuracy of the calculation is 
estimated to be of the order of 5%. This is mainly due to the strong sensitivity of the 
cross section to the nuclear model. The numerical accuracy of the calculations is of the 
order of 1%, and the cutoff dependence of the MEC is of the same order. 

8.4- Reactions with systems with A > 4 

Different from the lightest nuclear systems the photodisintegrations of nuclei with A > 4 
have not yet been studied in microscopic calculations with realistic nuclear forces, but 
LIT calculations with semirealistic NN interactions have been carried out |129[ll30lll31j . 
with the MT-I/III, MN, and AV4' NN potentials. Different from the two former 
potential models the AV4' interaction includes a nonzero force also for the relative P 
waves and also other odd waves. For the three- and four-nucleon systems such simple 
interactions lead to rather realistic total photoabsorption cross sections [951 1119j . thus 
also for A > 4 cases one can expect that they lead to similarly good results. 

In figure [28] we show the ^He and ^Li photoabsorption cross sections calculated with 
the above mentioned semirealistic potential models. Particularly interesting is the ^He 
case, which exhibits two separate cross section peaks. In a cluster model picture the 
low-energy peak can be explained by excitations of the surface neutrons (soft mode or 
'pigmy resonance') leading to a final state with an a-particle and two neutrons, while 
the second peak corresponds to the classical electric dipole resonance (relative motion 
of neutrons against protons) with a break-up of the a-core. A similar double-peak 
picture is not found for ^Li, even though a disintegration to a final state with an a- 
particle and a^So np pair is possible. However, such a transition is not energetically as 
well separated as the A^A^ + a channel and the a break-up channel in ^He. Note that a 
similar two-body break-up in case of ^He into the ^H-^H channel is not induced by the 
electric dipole operator. In figure [29] we show the comparison with experimental data. 
Compared with the other two interaction models the AV4' potential leads to a better 
agreement with experiment at lower energy both for ^He and ®Li. In fact for ®Li one 
finds a rather good description of the data up to 10 MeV, while at higher energies more 
precise data are needed in order to draw more definite conclusions. The agreement of 
theoretical and experimental low-energy results is less satisfying for ^He. On the other 
hand one has to take into account that the ^He photodisintegration cross section has not 
been determined from a direct measurement, but extracted from the Coulomb excitation 
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Figure 28. Total photoabsorption cross section of the six-body nuclei for various 
semirealistic NN potentials: AV4' (solid), MN (dashed), MT-I/III (dotted): (a) ''Li 
and (b) ^He. 



in the field of a lieavy nucleus. 

Now we turn to the ^Li photodisintegration. In figure [30] we show the cross section 
calculated with the AV4' potential in comparison with experimental data. One finds 
quite a good agreement in the whole considered energy range. The data exhibits a 
few fluctuations, but one would need more precise data in order to decide whether 
such structures in the cross section really exist. We would like to point out that 
the experimental ''Li total photoabsorption cross section had been determined via flux 
attenuation measurements in forward direction and not by summing up the cross sections 
of the various break-up channels. In the near future similar measurements for ^He and 
^'^Li are planned at MAXLab (Lund) and will hopefully also be performed at the HI7S 
facility (Duke). 

9. Summary 

We have reviewed the theory, tests and applications of the Lorentz integral transform 
method for perturbation induced reactions into the many-body continuum. The 
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Figure 29. As in figure [28] but experimental data are also shown: (a) ^Li{'j,n) data 
from |132j (solid squares) and ^Lz(7,^H)^He data from |133j (open circles) and |134| 
(solid triangles) added to the (7,n) data, (b) ^He data from |135j (soHd circles) 
(theoretical results convoluted with instrumental response function). 



method is a novel approach, which allows the microscopic calculation of inclusive 
and exclusive cross sections without using complicated many-body continuum wave 
functions. Since its proposal in 1994 the method has been applied to quite a number of 
electroweak processess with few-body nuclei, among them are the first calculations 
with a realistic nuclear interaction (NN and 3N forces) for the total ^H, ^He and 
^He photoabsorption cross sections [HSl I119j and the inelastic neutral current neutrino 
scattering off ^He [98] . Even reactions into the six- and seven-body continuum have been 
considered |129l 11301 1131] treating the many-body continuum interaction in a rigorous 
way. We would like to point out that at present there are no alternative approaches 
allowing similar ab initio calculations for A> 4. 

10. Acknowledgments 

W.L. and G.O. thanks the Department of Physics of the George Washington University 
for the warm hospitality. The work of N.B. was supported by the ISRAEL SCIENCE 



CONTENTS 



71 



2 



O 



1 


1 1 1 1 

AV4' — 
Ahrens et al ^ 














/ 1 


1 1 1 1 



o 



20 40 60 

CO [MeV] 



80 



100 



Figure 30. Total photoabsorption cross section of ^Li: tiieoretical result with AV4' 
potential (solid) and experimental data from [136j (circles). 



FOUNDATION (grant no. 361/05). V.D.E. acknowledges support from the RFBR, 
Grant No. 07-02-01222a, and the Russian Ministry of Education and Science, Grant 
No. NS-8756.2006.2. 



Appendix A. 

The correlated HH method ( CHH) 

In this appendix the CHH method is described. Not only the general idea of the 
method is summarized, but we also describe the details of its application, which has led 
to the results presented in section [81 

The basis functions we use are antisymmetric with respect to particle permutations. 
Such a basis is combined in a simple way from hyper spherical harmonics y belonging 
to a given type of permutational symmetry and from spin-isospin functions 6 belonging 
to a conjugated type of permutational symmetry. The set of 3^ differs from the 
unsymmetrized set (5.7) and we denote them y^pLMLa- They are components jj of 
a given irreducible representation [f] of the permutation group and they have a given 
grandangular momentum K, orbital momentum L with projection M^. The label a 
enumerates the various HH with the same K, L, Ml, [f] and fi. The 3^^^^^ q, may be 
constructed applying the Young operators to the simple of (5.7). This gives them 

in the form of linear combinations of the 3^[7^](Pfi), where P are particle permutations. 
These HH must be combined with spin-isospin functions OsmsTMt follows 



'KL,a 



(A.l) 



The spin-isospin functions have given spin S", Ms and isospin T, Mf quantum numbers 
and belong to the conjugated representation [/] of the permutation group. The spin- 



CONTENTS 



72 



isospin function with a /i value is conjugated to the spatial function with given /i so 
that the functions flA.l|) are antisymmetric with respect to particle permutations. They 



possess a given total momentum J and its projection M. When solving the dynamic 
equations, products of the functions (lA.ll) and hyperradial functions Rpi^p) can be taken 
as spatial basis functions x to expand |\E') in coordinate representation. For brevity we 
shall write the total basis functions as 

7 

where 9^ are the spin-isospin functions of the above mentioned form, and Xi-y ^"^^ the 
conjugated spatial functions. 

The short range repulsion at small distances, present in the central components of 
conventional NN interactions, slows down the convergence of the HH expansion (]A.2p 
considerably. To speed up the convergence one may include correlation factors in the 
basis functions. In [1^ it has been shown that, in the case of central NN potentials, the 
inclusion of a Jastrow correlation factor uj = Y[i<j fi'^^ij) the HH expansion accelerates 
the convergence drastically. To compensate for the strong short range repulsion potential 
that is present in the A-body dynamic equation the pair correlation function f{rij) 
should approximately be a solution of the two-body equation with this potential. When 
the mentioned simple Jastrow correlation factor is employed this may be fulfilled only 
for s-states of the NN pair and only when the triplet and singlet potentials V^i{rij) and 
Vi2,{rij) are close to each other in the region of short-distance repulsion. To go beyond 
these frames a more general ansatz has been proposed |137l [95] : 

|^)=a;^Q0i, (A.3) 

i 

where (pi are the HH-spin-isospin-hyperradial basis functions (1A.2I) and uj is the 
correlation operator of the form 

^ = n E fsT{r,,)PsT{t3)- (A.4) 

i<i S,T 

Here Psxiii) are the projection operators onto nucleon pair {ij) states of spin S and 
isospin T and S is the particle symmetrization operator. The inclusion of S is necessary 
since Psri'ij) for different pairs {ij) do not commute with each other. Operators of a 
similar type have been used in direct variational calculations of light nuclei, e.g. in |84j . 
In ( 1A.3|) |\E') stands for the initial state I^E'o)) or for the LIT function |\E'). 

This form of the correlation operator is easily incorporated into the calculations by 
first constructing the set of coefficients {9y\u}\9^) and then writing 

\^iI) = J2c^k (A.5) 

i 

where 

0i = ^Xi,7'^7'' Xi,7' = X](^7'l^l^7)Xi,7- (A.6) 
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The pair correlation functions entering (lA.4p are constructed as follows. At r < rg 
(r = Tij) they are solutions of the two-body equation 

\cPf 2df + 



M 



+ [V^Ar) + V^5tW]/ = (A.7) 



where / = fsr- Here I = for even states and I = 1 for odd states of a nucleon pair. 
The potential Vgj, is the central component of the NN interaction. The potential 
takes approximately into account the other components of the NN interaction and also 
the interaction produced by other nucleons. (In the case S = T = 1 the corresponding 
p-wave NN interaction is j-dependent which is disregarded in (]A.7p .) The potentials in 
(I A. 71) should ensure that /'(ro) = at some point tq. Equation (1A.7I) is solved in the 
range < r < tq and at r > tq one sets /(r) = /(ro) = 1. 

In the A = 3 calculations with realistic NN+3N forces presented in this review the 
procedure described above has been applied. For even NN states the potential Vg^ was 
disregarded in (]A.7p . In the A = 4 calculations with central NN forces the simple Jastrow 
correlation factor has been used with an average potential equal to l/2[V3i(r) + Vi3(r)]. 

In order to diminish the number of basis states in the correlated HH calculations 



one may optimize the potentials in ( ]A.7[) . In addition, one may apply different 
correlation operators ( ]A.4[) to different parts of |\E') or |\E'o)- The correlation functions 
may also be obtained from equations of (]A.7I) form for finite energies. These energies may 
be fixed from the condition that f^j. is continuous, i.e. fsTi^o) = [13], or in another 
way. In the ground state case the parameters determining fsrif^) niay be optimized 
from the minimum energy requirement or from the condition 

min {{{H - Eo)^o\{H - ^o)^o)}. (A.8) 

In the integral transform case these parameters may be chosen from the condition 

min {{{H - E)^ - Q\{H - E)^ - Q)} (A.9) 

for some typical E values, where \Q) indicates the source terms in the inhomogeneous 
equations (see section [2?T|) . In general, conditions of flA.8|) and flA.9l) form may be used 
to fix non-linear parameters in variational few-body calculations that are not based 
on the minimum energy principle. In particular in variational methods for calculating 
reactions the condition 

{{H - E)'^^\{H - E)"^^) = mm (A.IO) 

may be used for this purpose. 

Besides correlations caused by the strong short distance repulsion in NN potentials 
long range tensor correlations also require to reach high K values in the HH expansion. 
In the A = 4 CHH calculations presented here, as well as in some of the A = 3 ones a 
selection of correlated HH has been used. In general only HH that have small In and Kn-i 
values before applying the Young symmetrization operators are to be retained [Ml I138j . 
This prescription is justified for non-correlated HH and in the case of the ground state, 
but seems to work in our case as well. A more detailed selection within these frames 
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has been done (see also |126j .) For example, in the case of A = 3 bound state it has 
been found that for high K values only those HH are to be retained that in the notation 
of section 15.11 had quantum numbers I2 = 2 and /i = before applying the Young 
symmetrization operators. 

Convergent results have been obtained in the calculations with the CHH method 
(forjudging the accuracy reached see [96] in A = 3 calculations). Also in |139[ I140| 1141] 
accurate results have been obtained in the A = 4 case with realistic nuclear forces in 
a CHH basis for the ground state and low-energy scattering. There, different Jastrow 
type correlation factors for different parts of wave functions have been applied. It is 
ineteresting to notice that, due to flexibility of the basis, accurate results e.g. for the a- 
particle ground state with realistic NN+3N interactions, have been obtained retaining 
only low K values (up to = 8). 

A final comment regards the practical way to calculate matrix elements of the 
Hamiltonian between CHH states. In the A = 3 case an analytic integration over 
the Euler angles that determine the orientation of the system as a whole has been 
performed |142] . The remaining three-dimensional integration has been carried out 
with regular quadratures. In the A > 3 case the Monte Carlo integration has been 
applied. 



Appendix B. 

The effective interaction HH method (EIHH) 

In the effective interaction approach [143t 11441 \1AS\ I146t I147j the lowest eigenvalues 
of an A-body Hamiltonian 

H^^^ = H^ + V, (B.l) 

is treated in the following way. The Hilbert space of H^^^ is divided into a model space 
and a residual space, through the use of the eigenprojectors P and Q of Hq^ which 
satisfy the relations 

[ifo,P] = [^o,Q] =0; QHoP = PHoQ = 0; P + Q = l. (B.2) 

The Hamiltonian H^^'^ is then replaced by the effective model space Hamiltonian 

fflA]eff ^ pjj^p ^ Py[A]effp (g g^ 

that by construction has the same energy levels as the low lying states of H^^^. In 
general the effective interaction appearing in flB.3|) is an A-hodj interaction, and its 



construction is usually as difficult as finding the full-space solutions of the A-body 
problem. Therefore, one has to approximate V^^^^ff . However, one must build the 
approximate effective potential in such a way that it coincides with the bare one for 
P — > 1, so that an enlargement of the model space leads to a convergence of the 
eigenenergies to the true values. The EIHH method is developed along these lines. 
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Tp = -7^\. Tk{p) = —^ (B.6) 



In the HH formalism the internal A-body Hamiltonian H^^^ is written as 

m=f, + fK{p) + V^^\pA). (B.4) 

where 

A 

V^^\p,h) ^Y.v^, (B.5) 

denotes the bare two-body potential and 

1 . ^ . ^ 1 ^1-1 
-A,, Tk[p)-——, 

are the hyperradial and hypercentrifugal kinetic energies, respectively. In the previous 
equation Ap is the Laplace operator with respect to the hyperradial coordinate p, while 
Ka-i is the hyperspherical grand angular momentum (see equations fl5.5p - fl5.6l) ). The 
hyperradial kinetic energy Tp and the residual Hamiltonian 

f6^\p)^fK{p)+V^^\p,(l) (B.7) 

are often considered separately. The Hamiltonian 7i["^l(p), often used as a starting 
point in atomic and molecular calculations, is called adiabatic |148j as the hyperradial 
coordinate p is a slow coordinate with respect to the hyperangles. In the EIHH 
method |146l I147j the effective interaction is calculated starting with this adiabatic 
Hamiltonian. The unperturbed Hamiltonian Hq is chosen to be Tk{p) with the 
hyperspherical harmonics y[KA-i] eigenfunctions. The model space P is defined as 
the complete set of HH basis functions with generalized angular momentum quantum 
number K < Kp, and the Q-space as the complete set of HH basis functions with 
K > Kp. The states are denoted by {\p) ; p = 1,2, ...,np} for the P-space and 
{\q) ; q = rip+i, np+2, ...,nQ} for the Q-space. Of course, in principle one has uq — > oo. 
In actual calculations one considers finite Q-spaces. However, uq must be sufficiently 
large. 

For each value of the hyperradius p an effective adiabatic Hamiltonian is constructed 

?^[^i^^^(p, n) = PfK{p)P + P v^^'^'^ip, n) P . (B.8) 

However, as already pointed out, the effective potential would be a complicated A-body 
interaction, therefore V^^^^ff is approximated by a sum of quasi-two-body terms 

V[A]eff^ 

i<j 

Due to the use of antisymmetric wave functions one only needs to calculate the effective 
interaction operator relative to one pair, since 

iVl-^^ff)^{J2vQ'") = ^^^{vTd'-i,) ■ (B.IO) 

i<j 

The quasi-two-body effective potential v^a{a-i) determined as follows. First for each 
value p of the hyperradial coordinate one defines a quasi-two-body adiabatic Hamiltonian 
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containing the hypercentrifugal kinetic energy and the bare potential between the last 
two particles (the prefix quasi is justified, because this Hamiltonian depends on the 
A-body coordinate p) 

'H^^\p\Oa~i,Va-i) = fK{p)+VA,{A-i){V2psmOA^i-fjA^i) , (B.ll) 

The Hamiltonian of (IB. Ill) is then diagonalized on the A-body HH basis. Such a 
diagonalization is easily performed since p is only a parameter in Ti'^l, (there are no 
derivatives with respect to p) and for each value of p the Hamiltonian Ti^'^^p ; 9a-i, Va-i) 
depends only on three variables. This is just due to our choice of the A-{A — 1) 
pair in (IB.lOp . The obtained eigenstates are denoted by \j{p)) as they are continuous 
functions of p. 

One proceeds applying the Lee-Suzuki |149[ I150[ 1151] similarity transformation to 
'H^'^^{p ; 9a-i, Va-i) in order to get the corresponding hermitian effective Hamiltonian 

7YP]^^^(p ; Oa-i, fjA-i) = U\p) T^t'l (p ; ^a-i, r^^-i) U{p) , (B.12) 

where 

Uip) = {P + uj{p)) ^ . (B.13) 

.JP{l + u{pyCo{p))P 

The operator u}{p) is obtained using the following property |149j 

iuip) =Qiuip)P. (B.14) 

The matrix u;(p) is calculated for each value of p taking the np states |j(p)) with the 
lowest eigenvalues. Each of these states leads to the following system of {uq — np) 
equations 

{q\jip)) = J2^q\^^p^\p){p\jip))- (B.15) 

p 

The npijiQ — np) matrix elements (g|a)(p)|p) are obtained by solving the equation 
system ( IB. 151) . Once the effective quasi-two-body Hamiltonian iiP^^ff is constructed, 
the effective potential is obtained by a subtraction of the hypercentrifugal kinetic energy 



V 



'2]eff 
A,{A-1) 



ii^'^'"{p)-fK{p) . (B.16) 



Using this p-dependent effective potential and taking into account flB.Sp - QB.lOp one 
solves the A-body problem with the effective Hamiltonian 

Him^ = fp + 7^1^]^^^ = f, + fK + $^i)!f (B.17) 

in the P-space. One repeats the procedure enlarging the P-space up to a convergence 
of the low-lying energies of the A-hodj system. 

We would like to emphasize the following points: 

(i) it is evident that U{p) ^ 1 for P ^ 1 and thus v^'^^'^-ff converges to the bare va,{a-i)] 
therefore the energy spectrum converges to the exact one; 

(ii) in such a construction of an effective interaction the hyperradius is a parameter 
rather than a coordinate, and vf^^^^ is determined for various fixed p values; therefore. 
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Figure Bl. Binding energy (a) and root mean square radius (b) of ^He for the 
MTV potential, as a function of the maximal hyperangular momentum K„iax- The 
asymptotic value has been indicated by a dashed line. 



while being a two-body interaction, it depends on the whole A-hodj system via this 
collective coordinate; 

(iii) there is an additional dependence of v^'^^^f^ on the quantum number Ka-2 of the 
residual system (see equations (19) and (20) of [146] ): 

(iv) via the operator U{p) the effective potential contains information about a 
large part of the PQ-space interaction, hence the convergence to the exact eigenvalues 
of H^^'^ is accelerated with respect to the normal HH expansion; 

(v) application of the transformation operator U{p) to the hyperradial kinetic energy 
operator yields corrections to Tp that go beyond the adiabatic effective interaction |152j : 

(vi) the convergence of the EIHH method can be further improved through the use of a 
3-body effective interaction |152] . 

The effectiveness of the EIHH method can be best demonstrated comparing the 
convergence rate of a normal HH expansion with that of an EIHH one for some typical 
observable. As an example in figure Bl such a comparison is shown for the binding 
energy and root mean square matter radius of a four-particle system interacting through 
the MTV NN potential. From the figure it can be seen that, whereas the EIHH 
method gives almost converged results for Kmax = 4, the HH calculation with the 
bare interaction misses about 1.5 MeV for Kmnr = 20. 
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Appendix C. 

The multipole expansions of the response functions 

If the Hamiltonian is rotationally invariant it is useful to expand the states Q and Q' 
in (12.11) over states possessing given values J and M of the total angular momentum and 
its projection. Then the whole calculation may be done in separate subspaces of states 
belonging to given J and M. Furthermore, the calculations are M independent. Here 
we shall present the required relationships for the case of longitudinal and transverse 
electro disintegration response functions of a nucleus. 

We start from (I6.13p . (16.14^ . omitting the indexes 'a' and 

r-^(g,^) = ^y^ E ^c^/(vI/o(Jo,Mo)|/5t(q)|vI/XJ,M)) 

° Mo,J,M 

X (vl/^(J, M)|p(q)|vl/o(Jo, Mo)) 5 ^EjiJ) - Eo + ^ - uj^ (C.l) 
r^(?,^) = ^y^ E jtdf{MJo,Mo)\P{ci)\^fUM)) 

° Mo,J,M 

X (xl/^(J,M)|r(q)|vI/o(Jo,Mo))5(i?/(J)-Eo + ^-^) , (C.2) 

where the rotational quantum numbers J and M have been put in evidence. One then 
substitutes the expansions 

p{q) = A7:Y,VpUq)y;^{^), (C.3) 

jm 

A=el,mag jm 

into (IC.ip . (]C.2|) and apply the operation (47r)^^ J (iq to both sides of the arising 
expressions. This operation does not affect the responses while in their right-hand sides 
it leads to a simplification due to the ort honor mality property of spherical harmonics 
or vector spherical harmonics. This gives 



r^(g,u;) = -p— ^ it df {^o{Jo. Mo)\p]j^ f{J, M)) 

x(vl/XJ, M)|p,„|vl/o(Jo, Mo))5 (^Ef-Eo + ^-u^, (C.5) 

r^q.u^) = ^l^E E it df{^o{Jo.M,)\ffi\^>f{J,M)) 

° A MQ,J,MJ,m-^ 

X (M/^(J, M)\f^MJo, Mo))5 (^Ej-Eo + ^-S^. (C.6) 
Now we express these formulae in terms of the quantities 

Q'jM = (p. ® *o( Jo))^M . Qmi = [Tf ® ^o( Jo)) . (C.7) 
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They play the role of source terms in the LIT calculation. The notation of (1C.7P type 
means the Clebsch-Gordan coupling of the tensors of ranks j and Jq to the tensor of 
rank J. Using the orthogonality property of this transformation we finally obtain: 

r\q, uj) = 5^(2 J + l){r% (C.8) 

in 



where 



A=el,mag Jj 



,2 



xS{Ej-E, + ^-uj), (C.IO) 



{r^yf = ^ ci/(QSMl^/(^,M))(vl/,(J,M)|gSM) 



X 

,2 



x<5(i?^-i?, + -|_-^). (C.ll) 

The quantities on the right-hand sides of (IC.lOp . (IC.lip do not depend on M. In fact 
one may note that 

(v1/^(J,M)|QSm) = (2^+1)"'/'(^/('/)IIp.II^oW), (C.12) 
(vI/^(J,M)|QSm) = (2J+1)-^/^(v&/(J)||T;^||v&o(Jo)), (C.13) 

where the right-hand sides include the conventional reduced matrix elements. 

Similar to (13. ip we consider the partial transforms {L^Yj and {L^Yf of the 
contributions ( IC.lOp . (IC.lip . Taking into account the completeness of the set \1'/(J, M) 
in the subspace of states with given J and M one calculates these partial transforms as 

(L^)S = {{H -Eo~aR + iai)-^Q'j,,\{H -E^-aR + zaj)-'Q%), (C.14) 

= {{H -Eo-crn + iair^Q%,\{H -Eo-an + icyiY'QfM) -(CIS) 

Similar to (1C.12P and (IC.ISP the quantities (]C.14p and (IC.ISP are M-independent. The 
total transforms of the responses ( IC.SP and ( 1C.9P are 

L'{q, an) = ,p-- 5^(2 J + l){L%{q, an), (C.16) 

L\<i^^R) = i^p-7 E Y.^'^j+m^rXq.^R)- (ci?) 

A=el,mag J] 

To obtain the responses (IC.SP and (]C.9P from these transforms the integral equations of 
the form (13. ip are to be solved. Alternatively, one may invert the single terms separately 
and then sum up the corresponding partial contributions to the responses (IC.SP 
and (1091) . 
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The source terms (IC.7P possess definite parities. If V is tlie parity operator tlien 

vpjj" ={~iyp,m. (C.18) 
vffj^ ={-iyfjL (C.19) 
= {-iy+^f;^^\ (C.20) 

Tlierefore dynamic calculations proceed in subspaces of states with given J, M and given 
parities P = P{j)Po where Pq is the parity of the ground state and P{j) is (—1)-' or 

{-iy^\ 
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